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CHAPTER I 


INTRODUCTION 

The parabolic reflector antenna is the most useful and widely used 
antenna type for communications purposes at microwave frequencies. A 
parabolic reflector illuminated by a feed horn located at the focus of 
the parabola radiates a so-called "pencil beam" pattern which is a 
requirement for point-to-point communication. The other application of 
reflector antennas is in the field of radio astronomy. Normally, high 
gain and low noise temperature are required of the antenna for this 
application. The focus-fed parabolic reflector has limitations for 
large D/A since the spillover from the feed is usually pointed toward 
ground which is a high noise temperature source. Consequently, dual- 
reflector antennas derived from the Cassegrain telescope have been used 
for high gain and low noise temperature purposes. 

The parabolic reflector antenna with a feed at the focus preserves 
a fixed aperture distribution and thus, has a fixed radiation pattern 
when the feed pattern and the geometry of the reflector are given. The 
introduction of a second reflector in a dual-reflector system such as 
the Cassegrain antenna which consists of a hyperbolic subreflector and a 
parabolic main reflector or the Gregorian antenna which consists of an 
elliptic subreflector and a parabolic main reflector allows for more 
control over the aperture distribution of the main reflector. Although 
there is an extra degree of freedom provided by the addition of the 
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second reflector, the aperture distribution still possesses a certain 
amplitude taper which limits the performance of the antennas. 

Recent developments in the technology of space communications 
require that the antenna gain be optimized; or that the power radiation 
pattern of the antenna illuminates an irregularly shaped target area, 
the so called "contour beam" antenna. For an antenna to have optimal 
gain, the illumination on the main reflector must be such that the 
resulting aperture distribution makes maximum use of the entire 
reflector surface. For an antenna to have a power radiation pattern 
that follows a prescribed function in the contour beam application, the 
aperture distribution must be specially designed. None of these are 
provided by the conventional prime focus fed parabolic reflector or 
dual-reflector such as the Cassegrain or Gregorian antennas. 

Most recently, the reflector antenna is being used as a way to 
generate a plane wave which illuminates a target or antenna for 
measurement applications in indoor facilities. Electromagnetic 
measurements require that the antenna or scatterer under test be 
illuminated by a plane wave with a uniform amplitude distribution. In a 
compact range this is achieved approximately by a focused parabola. The 
non-uniform amplitude distribution caused by the taper of the feed-horn 
and the space attenuation limits the size of antenna or scatterer that 
can be measured. 

In the past 20 years, considerable research has been conducted in 
the area of dual-reflector antenna shaping to provide the capabilities 
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that the conventional reflector antennas are not able to achieve. These 


capabilities include: 

1. Maximize the illumination efficiency over the reflector 
aperture so that the antenna gain is optimized. 

2. Maximize the spillover efficiency so that most of the microwave 
power radiates in the desired direction and the power radiated 
to the other directions are minimized in order to reduce the 
interference between adjacent antennas and antenna noise 
temperature. 

3. Generate particular aperture distributions so that the 
radiation patterns of the antenna follow prescribed functions. 

4. Generate a plane wave which has uniform amplitude and phase 
distributions over most of the aperture of a compact range 
reflector so that the size of antenna and scatterer under test 
can be increased. 

A new approach of dual-reflector shaping is proposed and studied in 
this research. The main effort is limited to two-dimensional reflectors 
although the approach is demonstrated for three-dimensional shaping of a 
dual-reflector with circularly symmetric surfaces. This new approach is 
based on the geometrical properties of Cassegrain and Gregorian 
antennas. The shaped surface equations are formulated and solved 
numerically. 

Most reflector analyses and designs are based on the principles of 
geometrical optics. The method of physical optics is also widely used 
to calculate the scattered fields from reflector surfaces. These 
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techniques are briefly reviewed in Chapter II. Also, the geometrical 
properties of conventional Cassegrain and Gregorian reflectors are 
addressed in this chapter. 

Chapter III discusses the formulation for shaping of dual- ref lector 
antennas when the primary source and required aperture distribution are 
given. This is the basis of this research. A set of algebraic equations 
are obtained. These equations have to be solved numerically. A 
numerical approximation for solving these equations is presented in the 
same chapter. 

Several examples are presented in Chapter IV to validate the 
formulation and solution of the surface equations. A conventional 
Cassegrain or Gregorian reflector and feed pattern are assumed as a 
starting point, and then the reflectors are shaped to obtain a 
prescribed aperture distribution. The scattered fields of the 
subreflector and main reflector for both conventional and shaped dual- 
reflector antennas are calculated and presented for comparison. 

One problem associated with compact range reflectors is the edge 
diffraction from the sharp terminations of the main reflector surface. 
This diffracted field enters the area where the antenna or scatterer are 
under test and interferes with the plane wave reflected from the 
reflector surface. This interference causes variations on both the 
amplitude and phase of the plane wave illuminating the target under test 
and reduces the accuracy of the measurements. The edge diffracted 
fields can be reduced by adding rolled edges [1,2] to the main reflector 
surface. This modification to the edges of the shaped main reflector is 
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discussed in Chapter V. Examples are also shown to illustrate the 
improvement of adding rolled edges. The scattered fields from the 
reflectors are calculated by the method of physical optics which results 
in false end-point contributions to the scattered fields. The method of 
correcting the end-point contributions [3,4] is also reviewed in the 
same chapter. 

Chapter VI demonstrates the application of the shaping processes to 
three-dimensional circularly symmetric dual-reflectors that are 
illuminated by a primary source with a circularly symmetric power 
radiation pattern and thus the secondary aperture distribution of the 
reflector is also circularly symmetric. 

A summary and conclusion of the study are given in Chapter VII. 
Suggestions for further studies are also addressed. Three appendices 
are included. Appendix A presents a direct application of the 
principles of geometrical optics to dual-reflector shaping. Appendix B 
provides the details of solving the surface equations which are 
formulated in Chapter III. Finally, Appendix C shows the calculation of 
the geometrical optics reflected fields for the conventional Cassegrain 
and Gregorian reflector antennas. 
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CHAPTER II 


THEORETICAL TECHNIQUES 


This chapter briefly reviews the techniques which are used in the 
analysis and synthesis of shaped dual-reflector antennas. The designs 
of microwave reflector antennas are basically based on the principles of 
geometrical optics. As to the pattern analysis of reflector antennas, 
there are many methods available such as Geometrical Theory of 
Diffraction (GTD) and Physical Optics (PO). In this study, the method 
of Geometrical Optics (GO) is used to calculate the incident fields and 
the induced surface currents on the reflector surfaces, and then, the 
method of Physical Optics is used to calculate the scattered fields from 
the reflector. These methods are briefly described in this chapter. 

The method of synthesis of shaped dual-reflector antennas used in 
this research is based on the ray geometries of conventional Cassegrain 
and Gregorian dual-reflector antenna systems. The characteristics of 
the geometrical optics reflected rays of these antennas are reviewed. 

Throughout this research, the electromagnetic field is assumed time 
harmonic, and the e^ Wt time dependence is understood and suppressed. 

All the reflector surfaces are perfectly conducting and exist in free 
space. 
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2.1 Geometrical Optics 


The electromagnetic field associated with the propagation of 
visible light is characterized by a very high frequency of the field or 
very small wavelength. A good first order approximation to the 
propagation of electromagnetic waves at very high frequencies may be 
obtained by neglecting the finiteness of the wavelength. The branch of 
optics which is characterized by the neglect of the wavelength is known 
as Geometrical Optics. Mathematically, Geometrical Optics is the 
leading term in the asymptotic high frequency solution of Maxwell's 
equations and is subjected to certain geometrical laws. An excellent 
reference to geometrical optics is the work by Born and Wolf [5]. Some 
basic properties of the geometrical optics field are 

1. The time-averaged electric and magnetic energy densities are 
equal, and they are each equal to one half of the total stored 
energy density. 

2. The average Poynting vector is in the direction of the normal 
to the geometrical wavefront, and its magnitude is equal to the 
product of the average energy density and the velocity v=c/n 
where c is the speed of light in free space, and n is the 
refractive index of the medium. In free space, n=l so that 
v=c. 

3. Geometrical optics is a ray tracing technique in which 
geometrical rays are the family of directed curves 
perpendicular to the geometrical wavefront. Geometrical rays 
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also represent the direction of energy flow at every point. In 
a homogeneous medium, the geometrical rays are straight lines. 

4. The geometrical optics field is a transverse electromagnetic 
field (TEM), i.e., the electric and magnetic field vectors at 
each point are orthogonal to each other and to the ray. 

5. In a homogeneous medium the polarization of the geometrical 
optics field remains constant along a ray. 

6. A ray optical field can be described in terms of an astigmatic 
ray tube as shown in Figure 2.1. The electric field at s 
associated with ray tube is given by 

l(s) = £(o) A(p lf p 2 ,s)e (2.1) 

in which t(o) is the field at the reference 0, A(p^,p 2 ,s) is the spread 
factor, and p^,p 2 are the radii of curvature of the geometrical wave 
front. The spread factor is given by 


A(p 1 ,p 2 ,s) = 

which results from the conservation of energy along the ray. Equation 
(2.1) expresses the intensity law of geometrical optics. Since the 
geometrical optics field is TEM, the magnetic field must also follow 
these laws. If one of the radii of curvature is infinity as in the two- 
dimensional case, the spread factor reduces to 


P-, P, 


12 


(p.+s)(p 2 +s) 


( 2 . 2 ) 
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Figure 2.2. Geometry associated with the reflection by a surface S. 


\ 



Figure 2.3. Induced surface current in the physical optics 
approximation. 
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r r 

s = distance from reflection point to observation point and P^»p 2 
are the principal radii of curvature of the reflected 
wavefront at the point of reflection Q R . The unit vectors 
associated with the dyadic reflection coefficient are defined 


ej|= unit incident vector parallel to the plane of incidence and 
perpendicular to the incident ray 
eM= unit reflected vector parallel to the plane of incidence and 


perpendicular to the reflected ray 


'i 


unit vector perpendicular to the plane of incidence where the 

plane of incidence is defined as the plane formed by the unit 
*i 

vectors s and n as shown in Figure 2.2. 


it r 

The general expressions for Pj and p 2 can be found in [6]. 
However, for some special surface geometries, p^ and p 2 can be obtained 
numerically or directly from the properties of the surfaces. When one 
of the radii of curvature goes to infinity, as in the two-dimensional 


case, the reflected electric field becomes 


2 r (s) = f i (Q R ) 



(2.5) 


For a two-dimensional surface which is infinitely long along the z-axis, 
the reflected electric field for a TE polarized incident wave (generated 
by a z-directed electric line source) is given by 
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f r (s) = E r (s)z = -zE 1 (Q r ) 


( 2 . 6 ) 


k« 

and the reflected magnetic field for a TM polarized incident wave 
(generated by a z-directed magnetic line source) is given by 


i? r (s) = H r (s)z = zH i (Q R ) 

One final note for the geometrical optics reflected field is that the 

law of reflection is satisfied; i.e., the angle of incident 6* and the 

angle of reflection 6 r are equal, and the incident vector s 1 , the 
A r 

reflected vector s and the surface normal n are co-planar. 

2.3 Physical Optics 

Physical optics (PO) is also a high frequency method which uses the 
integral equation representation, along with the physically reasonable 
high frequency assumption that the scattered field from one point on a 
body to any other point is insignificant compared to the incident field 
strength. This method is based on the approximation of Stratton-Chu 
Equations for the scattered field [7] and states that for scattering 
from a perfectly conducting body in free space, the scattered fields are 
given by 


-jks 


r 

p.+s 


(2.7) 
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r ■ - Jk o z o Js p s * + ri < J s • V H ds 

K. 


( 2 . 8 ) 


and 


& = J g x 7 *) ds. (2.9) 

Note that <J> is the free space Green's function, and 3 is the induced 

s 

surface current on the scatterer surface as shown in Figure 2.3. In 
general, the induced surface current 3 is given by 

3 « n x j| T0TAL « n x (if 1 + II s ) (2.10) 

in which j| T0TAL i s the total magnetic field on the surface of the 
scatterer, if 1 is the incident magnetic field on the surface, and n is 
the surface normal. Substituting Equation (2.10) into Equations (2.8) 
and (2.9), it is apparent that the unknown field if 8 appears on both 
sides of the equation so that a coupled set of integral equations is 
obtained. Although one can solve if s by the Method of Moments, it is 
very tedious and inefficient since the size of the scatterers is usually 
very large in terms of a wavelength in the high frequency region. 
Consequently, from the boundary conditions at the perfectly conducting 
surface, one obtains that 

n x t 1 = n x I s on surface S. 

The induced surface currents are approximated by 
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'2 n x if* on the lit side 
,0 on the shadow side 


( 2 . 11 ) 


and the scattered fields 2 s and if s can be calculated by Equations (2.8) 
and (2.9). This is the so-called "physical optics" approximation. The 
approximation to the induced surface currents are valid when the 
transverse dimensions of the surface, the radii of curvature of the 
surface, and the radii of curvature of the incident wave front are all 
much larger than a wavelength. 

For the three-dimensional case, the free space Green's function is 

-jkr 

where r is the distance from the source point to the field point 

4 nr 

and when r»X, Equations (2.8) and (2.9) become 


f S 


and 


jk Z 
_ o o 

4n 








ds. 


( 2 . 12 ) 


(2.13) 


For the two-dimensional case, the free space Green's function is 
ij H^^(kp) where H^^(kp) i s the zero** 1 order Hankel function of the 

2 n< * kind, and the surface integral in Equations (2.8) and (2.9) should 
be replaced by a line integral. There are two polarizations in the two- 
dimensional problem, one is the Transverse Electric (TE) case, and the 
other is the Transverse Magnetic (TM) case. In the TE case, the 
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incident fields are generated by an electric line source along the z- 
axis; while in the TM case, the incident field is generated by a 
magnetic line source along the z-axis. 

For an electric line source, the incident electric field is i^-zE* 
and the incident magnetic field is as shown in Figure 2.4. The 

physical optics induced surface current l g is given by 

^=2nxit i *2nx #1^ = z CH^ » zJ . (2.14) 

s ♦ ♦ z 

Thus, from Equation (2.8) with the surface integral replaced by a line 
integral, one finds that 


I s - 


E S z * 


-jk Z f z J * + K (J z • 7) 7* dl 
O O J Z z 

over 

illuminated region 


(2.15) 


1 * 

But As (J z • 7)7i|> = 0 since there is no variation along z-axis, i.e., 
k Z z 

h * 0> thus 


K ■ ~K z o R • SJ dl - - t 2 I J z h o 2> < k ’> dl 


(2.16) 
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For a magnetic line source, the incident magnetic field is H =z Hz; 
thus, the induced surface current on the lit side is given as 

5 * 2(n x it*) 

s 

3 = 2(n x z Hz) = 1 C Hz 
s 

or 


J 


s 



(2.17) 


where 1 is the unit tangent vector along the surface of the scatterer as 
shown in Figure 2.5. The scattered magnetic field is then given by 


3 s - Jl U 1 x V* dl 


(2.18) 



+ (-H< 2) (k P )) 

V* = - P H< 2) (kp) 


18 


MAGNETIC 
LINE SOURCE 


Figure 2.5. 


TM incidence in two-dimensional case. 



Consequently, one obtains 


if 8 = - J x J x i x (+ £y p h[ 2) (1cp)) di 

or 

S® - - £y J (i X *p) J x H< 2) (kp)dl (2.19) 

and the scattered magnetic field is given by 

if 3 = zH 8 = - z | z • (i x p) H^ 2) (kp) dl. (2.20) 

When the observation point is in the far-zone of the induced current the 
large argument approximation of the Hankel function can be used and 
Equations (2.16) and (2.20) become 


E 


s 

z 



e jn/4 



dl 


( 2 . 21 ) 


H 


s 

z 



z*(l x p) ^ 

.lit 



dl 


( 2 . 22 ) 


These two equations are used in this research to calculate the scattered 
fields from the shaped reflector antennas. 


20 


2.4 Conventional Cassegrain Reflector Antenna 


The conventional Cassegrain reflector antenna is a two-reflector 
antenna which is designed in the form of a hyperbolic subreflector and a 
parabolic main reflector based on the principle of the Cassegrain 
optical telescope [8]. The geometry of the Cassegrain system is shown in 
Figure 2.6. The feed of the antenna is usually located at the real 
focal point of the hyperbolic surface while the focus of the main 
reflector is matched with the virtual focal point of the hyperbolic 
surface in a focused system. The surface of the subreflector can be 
described by 



and the surface of the parabolic main reflector by 


x = ¥=- - (F -F ). 

4F v m c' 
m 


(2.24) 


The various parameters associated with the reflector surface are defined 
below: 

F c = distance between the two focii of the hyperbolic surface 
F = distance between the real focus and the vertex of the 
hyperbolic reflector 
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L = F -F = distance between the virtual focus and the vertex of the 
v c 

hyperbolic surface 


e = eccentricity of the hyperbolic surface = p — — 

c“ v 

2a = length of the transverse axis of the hyperbolic surface = F/e 

2b = length of the conjugate axis of the hyperbolic surface 

12 / 

= a e -1 


Pattern analyses of the Cassegrain antenna have been well 
established such as the one in 19]. One can refer to this reference for 
more detail. However, the behavior of the reflected rays from the 
subreflector and main reflector will be summarized here. For a focused 
hyperbolic surface illuminated by a source located at the real focus, 
the reflected rays from the surface appear to emanate from the virtual 
focus. Thus, the reflected field caustic distance is equal to the 
distance between the virtual focus and the point of reflection on the 
subreflector. The reflected rays from the main reflector are then 
parallel to the reflector axis since the incident rays radiate from the 
virtual focus of the subreflector; i.e., the focus of the parabola. The 
ray geometry is given in Figure 2.7. This phenomenon is the basis for 
the dual-reflector shaping of this research. 
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I 

I 

Figure 2.7. Ray geometry of the Cassegrain antenna system. 
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2.5 Conventional Gregorian Reflector Antenna 

Similar to the Cassegrain reflector antenna, the Gregorian antenna 
is also a two-reflector system but, instead of the hyperbolic 
subreflector, the subreflector is an elliptical two-dimensional surface 
and the focus of the main reflector lies between the two reflectors as 
shown in Figure 2.8. The elliptic surface has two focii and the focus 
of the main reflector is usually matched with focal point #2 while the 
primary feed is located at focal point #1 in the focused case. The 
surface of the elliptic subreflector can be described by 


x = a 



(2.25) 


and the equation for the surface of the main reflector is the same as 
Equation (2.24). The parameters of the antenna are 


F = distance between the two focii of the elliptic surface 
c 

F = distance between focal point #1 and the vertex of the elliptic 

reflector 

L = F -F 
v c 

F 

c 

e = eccentricity = f _2 l — 

c v 

2a = length of the major axis of the elliptic surface = F c /e 

I 2 -1 

2b = length of the minor axis of the elliptic surface = a U-e . 
Note that F <F so that L is negative and e<l. 

C V 
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The ray geometry of the Gregorian antenna is shown in Figure 2.9. 
When illuminated by a source located at focal point #1, the reflected 
rays from the elliptic subreflector pass through focal point #2 and the 
reflected rays from the main reflector are again parallel to the 
reflector axis since the incident rays emanate from focal point #2. The 
caustic distance of the reflected ray from the subreflector is the 
distance between focal point #2 and the point of reflection on the 
subreflector surface. However, this caustic distance is negative since 
the reflected rays pass through it. 
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CHAPTER III 


SYNTHESIS OF TWO-DIMENSIONAL DUAL-REFLECTOR ANTENNAS 


3.1 Introduction 

The design of dual-ref lector antennas, including the conventional 
Cassegrain parabola-hyperbola and Gregorian parabola-ellipse antennas, 
are all based on the principles of geometrical optics. Since 
geometrical optics is a high frequency method, the reflectors must be 

t 

large and have a large radius of curvature compared to the wavelength. 
The principles of geometrical optics state that: 

1. Snell's law [10] must be satisfied at each reflector, that is, 
the incident ray, the reflected ray, and the surface normal at 
the reflection point on each reflector must be coplanar, and 
the angles of incidence and reflection are equal. 

2. Power flow along each differential tube of rays remains 
constant, even when the tube undergoes a number of reflections 
(conservation of power [11]). 

3. Ray directions are normal to the surface of constant phase, and 
this condition is maintained after a number of reflections 
(theorem of Malus [12]). 
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Conventional Cassegrain and Gregorian reflector antennas are two 
examples which satisfy the above principles. The geometries of these 
two antennas are shown in Figures 2.6 and 2.8 of the previous chapter. 

As discussed previously, with the source which has a spherical wavefront 
(or cylindrical wavefront for the two-dimensional case) and is located 
at focus #1, the reflected rays from the subreflector which satisfy 
Snell's law of reflection appear to emanate from the second focus and 
also have a spherical wavefront (or a cylindrical wavefront). Since 
the second focus is also the focus of the parabolic main reflector, the 
reflected rays from the parabola which also satisfy law of reflection 
form a uniform plane wave and are parallel to the axis of parabola. It 
can be shown that the power contained in the ray tubes before and after 
reflection are equal in order to satisfy the law of power conservation. 
Also, the incident and reflected rays from each reflector are normal to 
their corresponding wavefront. 

In general, the geometrical optics designs of dual-reflector 
antennas must satisfy the above principles either mathematically or 
physically. Galindo [13] and Kinber [14] used these principles to set up 
two first-order ordinary differential equations which can be solved 
exactly, if possible, or numerically to obtain a couple of circularly 
symmetric reflectors for realizing a prespecified aperture amplitude 
distribution. Green [15] and Williams [16] also used these principles 
to design the subreflector and then correct the main reflector to 
achieve the desired uniform phase criterion at the aperture plane of 
main reflector. A direct application of principles of geometrical 
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optics in dual-reflector antenna shaping is given by [17] and is 

discussed in Appendix A. For offset dual-reflector antennas, these 

principles are again used to obtain partial differential equations which 

are then solved numerically to obtain the surfaces [18,19,20]. Westcott 

* 

[21] used complex coordinates to set up a Monge-Ampere type partial 
differential equation, and solved this equation as a boundary value 
problem to synthesize a dual-reflector antenna. He also used the same 
approach to synthesize a single reflector for relating a given feed 
pattern to a prescribed far-zone pattern. In the following section, a 
new approach which is based on the knowledge of ray behavior of the 
conventional Cassegrain and Gregorian antennas is presented. Instead of 
setting up differential equations or heuristically designing the 
subreflectors, a set of nonlinear algebraic equations are obtained and 
solved numerically to obtain the surfaces. Also the caustics of the 
reflected rays from the subreflector can be obtained and examined. 

3.2 Method of Synthesis 

As discussed in the previous section, a focus-fed Cassegrain or 
Gregorian reflector antenna will generate a uniform plane wave on the 
aperture plane of the main reflector. However, for a given primary 
source, the amplitude distribution of the aperture field preserves a 
certain amplitude taper which is fixed when the source and the 
parameters of the reflector surfaces are specified. From geometrical 
optics, it is known that the amplitude of the reflected field is 
dependent on the space spread factor or the location of the caustic of 
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the reflected wavefront. In order to realize a given aperture amplitude 
distribution for a given primary source, the incident rays on the main 
reflector which are also the reflected rays from the subreflector have 
to be different from those for the conventional Cassegrain and Gregorian 
reflectors. Also, in order to generate a uniform plane wave, the main 
reflector locally must have the properties of a pure parabola with its 
axis parallel to the ray direction of the desired plane wave. Since the 
caustic of the reflected ray from the subreflector is also the focus of 
the main reflector, this caustic will not be a fixed point as in the 
conventional Cassegrain and Gregorian antennas so that the amplitude 
distribution of the aperture field can be controlled by moving the 
caustic for each ray reflected from the subreflector. The formulation of 
this synthesis method will be discussed in detail in this section for 
two-dimensional reflectors. The general geometries are given in Figure 
3.1 for both concave and convex subreflectors. The coordinate origin is 
assumed at the point where the primary source is located. Although these 
formulations are basically for two-dimensional geometries, they can be 
easily modified for shaping of three-dimensional dual-reflector antennas 
with circular symmetry. 

The basic idea is that for a differential tube of rays emanating 
from the primary source, the subreflector section between the 
intersection of the incident ray tube and the subreflector (i.e., the 
surface between ( x 10 >y 10 ) and (x^y^) is assumed to be either a 
hyperbola or an ellipse with the two focii at the coordinate origin and 
the caustic point ( x c »y c )* T ^ e corresponding section on the main 
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reflector (i.e, the surface between ( x 2o»y2(P an< * ^ x 2’^2^ * s assume< * t0 

be a perfect parabola with focus at (x c »y c ). Also shown in Figure 3.1 

is (x ,y ) which is the caustic of the previous reflected ray tube, 
co co 

Based on these assumptions, one can set up a set of equations and solve 
these equations to obtain a solution of (x c ,y c ), (x^,y^) and (Xg^)* 
These solutions will be accurate as long as the width of the incident 
ray tube from the primary source is small. The formulations for the 
surface are given below. Note that the "old" ray in Figure 3.1 is 
assumed to be known; i.e., (^O^IO^’ ^ x 20’^20^ an< * ^ X co’^co^ are ^ nown * 
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(a) Convex subreflector 
Figure 3.1. Geometry of shaped dual-reflector. 
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(b) Concave subreflector 


Figure 3.1. Continued. 
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3.2.1. Formulation 
A. Subreflector 

The canonic form of a hyperbola as shown in Figure 3.2 can be 
written as 



where 2a is the transverse axial length, 2b is the conjugate axial 
length, and b is related to a by 



Note that F is the distance between the two focii F' and F, and 
c 

2a * FP - FP . 

The canonic form of an ellipse as shown in Figure 3.2 can be 
written as 



where 2a is the major axial length, 2b is the minor axial length, 
axial length can be related by 


The 
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2 2 
4a -F 

c 


JL , 2 
F -4a 
c 


Also, one finds that 


2a « F'P + FP 

o 

By observing that b in both the concave and convex cases has the same 
expression except for the sign, one can write, in general, the surface 
as 



2 

where b <0 for an ellipse 
2 

b >0 for a hyperbola 

2a = F'P + FP for an ellipse, and 

2a = F'P - FP for a hyperbola 

From these canonic forms for the hyperbola and the ellipse, the 
shaped subreflector section between ( x iQ' y io^ an< * ( x i» y i) which is 
assumed to be a hyperbola or an ellipse can be formulated by proper 
rotation and translation of the corresponding canonic form. As shown in 
Figure 3.3, the distance F £ between the two focii for this section is 


38 





given by 


„ / 2 2 . 1/2 

F c - < x c + V ■ 


The transverse axial length is given by 


/ 



v 



+: ellipse 
hyperbola 


where ( x io ,y l(p * s assumel * t0 known a priori. The subreflector 
surface section is then defined by 



but x" = x'cos© + y'sin© and y" = -x'sinO + y'cos0. 


Then, x' and y' are related to x and y as 


x' 



and 


y' 



The surface in the x-y coordinate is given by 
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^2 j( x “ 2^) cos0 + (y ~ 2 ^) sine • 


2 < - 2 ^J sin6 + |y - cosej = 1 


which can be further simplified to 


1 ( F d 2 4 { ) 2 

~2 jx cos6 + ysin6 - l 2 2 — | x s * n ® ~ ycosel = 1 


F - 4a 
c 


since F cose = x and F sin6 = y . 
c c c ■'c 


The unknown subreflector surface point (x^,y^) must thus satisfy the 
following equation: 


1 / F y 

~2 tacose + yj^sine - ^1 


2 2 
F -4a Z 
c 


(x^sine - y^cosej' 


(3.1) 


where 


„2 ( 2 2^ 

F c ■ K + y cJ 


4a J - 


( r iw \ 2 

(2 2 W ( )2 ( )2 

- ' l X 10 + y loJ ± l x 10- x cj + l y l(T y cJ ’ 

^ L J / 


+: ellipse 

hyperbola 
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B. Main Reflector 

The canonic form of a parabola, as shown in Figure 3.4 can be 
described as 



where F is the focal length of the parabola and can be expressed as 
F = j (1 - cos a)p . 

Now, by assuming that (x 2 o ,y 20* and ( x 2’ y 2 * form a P arabola with the 
focus at (* c ,y c ) as shown in Figure 3.5, the following expression is 
obtained to describe the main reflector surface: 

* - M -?F (y - y J 2 

where 

F = | (1 - cos a) p 
or 

\ ^~ sec a o K V X 20 ) 

F = 1 l -cos a 0 

T2 sina ^ y c -y 20^ 

0 
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y 2 

PARABOLA : X = ~T=- 


Figure 3.4. Canonic form of parabola. 
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Therefore, the new main reflector surface point (x 2 ,y 2 ) must satisfy 

* 2 - <V P) ” h <y 2- y c )2 - (3 ' 2) 

Equations (3.1) and (3.2) give two equations with 6 unknowns 
( x l , yi ,x 2 , y2 ,x c’^c^ ’ ^ our more e< I uat i ons ar ® required to solve for all 
the unknowns. Basically equations (3.1) and (3.2) satisfy the theorem 
of Malus which has been discussed earlier. The other two conditions of 
geometrical optics will be used to obtain four more equations. 

C. Snell's Lav of Reflection 

At the point of reflection, the angles of incidence and reflection 
must be equal, and the incident ray, reflected ray, and surface normal 
must be coplanar. Since the antenna is assumed to be in a homogeneous 
medium, the incident and the reflected rays must be straight lines as 
shown in Figure 3.6. Consequently, one finds that 

y 1 = (tan$)Xj = m ^ (3.3) 

and 

(x 2 -x 1 )(y c -y 1 ) = (y 2 -y 1 )(x c -x 1 ) (3.4) 

The angles of incidence and reflection will automatically be equal since 
the reflection point (Xj,y^) is on either the hyperbolic or the elliptic 
subreflector. Snell's law of reflection will also be satisfied 
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Figure 3.6. Geometry for describing the relations of various points. 
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automatically on the main reflector. One more equation can be obtained 


from Figure 3.6 by considering that the intersection point of the "old" 
ray with the "new" ray must be ( x c »y c )» 

x -x x on -x in . 

— — — , JO-JO . (3.5) 

V y co y 20- y 10 tan “o 

which can be rewritten as 


and 


x + Ap cosa 
co o 


y 


c 


y co + Ap sin« o 


Note that Ap is the unknown distance between (x ,y ) and (x ,y ) and 

c c co co 

can be positive or negative. 


D. Conservation of Pover 

The major purpose of shaping the reflector is to redistribute a 

given primary source pattern into the prescribed secondary aperture 

distribution. Thus, conservation of power has to be satisfied. 

Consider a primary source with a power density pattern F($) and with its 

axis tilted with respect to the positive x-axis. Note that the tilt 

angle is zero for a symmetric reflector. The power between 

and $=$ of this primary source is to be distributed over a secondary 
max 

aperture between y=y m j n and y = y max with an aperture distribution CI(y). 

The angles 4, 4> . and $ are also measured with respect to the 
6 T ’ min max r 
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positive x-axis. The constant C is determined from the conservation of 
power which states that 


^max 

F(40d4> 

, ^min 


v 

•'max 

C I(y)dy 

y . 

J'min 


max 

6 

J min 

F(*)d* 



•'max 


y t 

i(y)dy 

J^min 



Note that in this study, the aperture blockage effects from the 
subreflector are not considered. It is assumed that the shaping process 
begins at the point where the incident ray from the primary source is 
along the feed axis and the surfaces above and below this 

point are shaped separately. Thus, the conservation of power relates 
the feed angle $ and the aperture point y as 

V 

C I(y)dy (3.6) 

y o 

where y Q is the point at the aperture which corresponds to the ray which 
emanates along the feed axis. The geometry of describing conservation 


F(*)d<> = 
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of power is shown in Figure 3.7. For a given y (or $), the 
corresponding $ (or y) can be solved from Equation (3.6). 

The equations obtained thus far are summarized as follows: 


^2 {xjCOsG + y^sinG j |x^sinG - y^cosoj = 1 (3.1) 


F -4a 
c 


x 2 - <V F) ' SF <y 2- y c )2 


(3.2) 


y 1 - (tan<|>)x 1 


(3.3) 


(x 2 -x 1 )(y c - y i ) - < x c - x i )(y 2" y i ) 


(3.4) 


x c x 20 _ x 20 x 10 = 1 

y c -y 20 ’ y 20 -y 10 = ,ami ( 

and 


(3.5) 


F(+)d+ 


max 


min 


F(*)d* 


(V, 


i(y)dy 


max 


I(y)dy 


min 


in which 


(3.6) 
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Figure 3. 



y =y, 


max 


y = y 0 


- X 


y =y, 


mm 


. Geometry of describing conservation of power. 
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{< ♦ y2 r 


(3.7) 


F » 



and 


+ : ellipse 

_J»yperbola 


F = 


2 ( 1 - sec “o) (V x 20) 


- 1-cosa 

1 o 

2 sina 


( y c- y 20) 


, if x 

» if y c * y 


20 

20 * 


Although there are 6 equations with 7 unknowns (Xj,yj,X 2 »y 2 »x c ,y c , <j>) 
one of these unknowns will be used as an independent variable so that 
the other 6 unknowns can be solved by the 6 nonlinear algebraic 
equations. 

It is also noted that the method of shaping discussed above will 
generate uniform phase on the aperture when the primary source has a 
cylindrical wavefront since the path length from the primary source to 
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the aperture plane will be the same for two adjacent rays which are from 
the same section of a Cassegrain or Gregorian surface. 

3.2.2. Solution for the Surface Equations 

It is obvious that to solve the equations obtained in the previous 
section will not be easy even by numerical methods since they are 
nonlinear. In this section, an approach is proposed to solve these 
unknowns . 

As mentioned earlier there are 6 equations with 7 unknowns? 
therefore, one of the unknowns must be chosen as an independent 
variable. The obvious choice will be either y ^ or f since once one of 
these two variables is specified, the other one can be obtained from 
Equation (3.6). This reduces the set of equations to five equations 
with five unknowns. It is also noted that the above equations involve 
"old" and "new" rays and in the very first stage, the "old" ray is not 
defined unless one provides a set of initial points {( x i 0 ’^lO^ » ( x 20 '^ 20 ^ 
and (x ,y )l. Once these initial points are set, the method which 

CO CO ) 

will be discussed next can be used to solve for x ^ , yi >x 2 ,x c , ^c an< * 
these newly generated points will be used as the "old" ray to solve the 
next "new" ray. This process is continued until the complete surfaces 
are generated. 

The procedure of solving the above equations is discussed as 
follows: 

1) From Equation (3.3), 


52 


(tan^Xj = 


(3.7) 


2) Substituting (3.7) into (3.4), one obtains 


x c y 2 - x 2 y c 
1 “ y 2- y c" m o (x 2- x c ) 


x, = 


(3.8) 


3) Substituting Equation (3.2) into Equation (3.8), one obtains 


4F X c y 2~ y c [ (y 2~ y c )2 + 4Fx c ~ 4f2 ] 

4F(y 2 -y c ) - %{[(y c -y c ) 2 + 4Fx c - 4p2 ] - 4Fx J 


(3.9) 


4) Substituting Equations (3.7) and (3.9) into Equation (3.1), an 
equation which involves x c »y c * m 0 » y£ obtained. The resulting 
equation is 


(vvc) 2 - 4a2 1 1 * "o) 


“cV'tll'rtF * 4Fx c - 4p2|k 


/ 

- (vVc) ( F c- 4a2 ) ' 4Fx c y 2- y c [( y 2 - y c) 


i2 +4FX - 4F 2 
c 


4P ( y 2-Wc) - m o[( y 2- y c) 2t4Fx c- 4p2 ]} 


' - l ( F c- 4a2 ) 2 ' 4F ( y 2- y c*Vc) -m o [( y 2 -y c) 2+4Fx c -4p2 ] 


(3.10) 
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Although Equation (3.10) has and F c involved, it is noted that 

both m Q and are known at this point and that F and F c are functions 

of x and y only, 
c c 

5) There are two equations, Equations (3.5) and (3.10), with two 
unknowns x c and y £ which can be solved. However, Equation (3.10) is 
again a nonlinear equation so that it is still very difficult to solve 
for x £ and y £ exactly. Consequently, a numerical method has to be used 
to solve for x £ and y £ . One can directly solve these equations 
numerically, or by observing that Equation (3.5) can be written as 

x = x + Ap cosa 
c co o 

and 

y ■ y + Ap sina (3.11) 

7 c •'co o 

in which Ap is the unknown distance between ( x co >y co ) an d ( x c >y c )* 

After substituting Equation (3.11) into Equation (3.10), only one 

\ 

unknown Ap is involved. But this equation is still not easy to solve. 

However, it is expected that Ap will be very small when the angle $ (or 

the aperture point y 2 ) is changed gradually (i.e., the increment of 4> is 

2 

very small) one can neglect the higher order terms of Ap such as (Ap) , 

3 

(Ap) , etc., in the newly obtained equation. After Ap is solved, all 
the unknowns can be solved. The detailed procedure of solving these 
unknowns is given in Appendix B. 
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CHAPTER IV 


EXAMPLES FOR SHAPING OF TWO-DIMENSIONAL DUAL-REFLECTOR ANTENNA 


The method of synthesis for dual-reflector antennas discussed in 
the previous chapter is actually quite general. For a given feed 
pattern, feed location, angular range of feed pattern which will 
illuminate the subreflector, and the desired aperture size and 
distribution, the reflector surfaces can be generated numerically 
provided that a set of appropriate initial points for the caustic, 
subreflector and main reflector is given. This set of initial points 
must satisfy a ray condition that the resulting aperture field obtained 
by connecting the rays through these points is equal to the desired 
aperture field. Two of the three initial points are chosen, the third 
one can be calculated. Generally, the initial points on the 
subreflector and main reflector are chosen and the initial point for the 
caustic is then calculated by the formula discussed in Appendix C. 
However, as it turns out, if the three initial points do not satisfy the 
above ray condition, the solution for the surface itself will try to 
correct the set of inappropriate initial points. This phenomenon can be 
seen from the examples given later in this chapter. 

Several examples of two-dimensional dual-reflector shaping are 
presented in this chapter to illustrate the method of synthesis 
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discussed in the previous chapter. In these examples a conventional 
Gregorian or Cassegrain antenna is used as a starting point for the 
reflector surface; and a magnetic line source with power radiation 
intensity (cos^) is used as the primary source. The goal of shaping is 
to modify the given reflector surfaces so that a prescribed aperture 
distribution with uniform phase is achieved at the aperture of the main 
reflector. The physical optics approximation is then used to calculate 
the scattered field from the shaped surfaces although the geometrical 
optics is used to compute the induced surface currents on the reflector 
surfaces. 

In each example, the location of a "central ray" at the aperture is 
assumed known where the central ray is the reflected ray from main 
reflector which corresponds to the incident ray coming along the axis of 
primary feed. This central ray is used to determine the initial points 
of shaping. The reflector surfaces above and below the central ray are 
then shaped separately. The initial points on the reflector surface can 
be determined by the intersection of the central ray with the initial 
Cassegrain or Gregorian reflector surface; and the initial point for the 
caustic is determined by the condition given in Appendix C. 

In addition, as mentioned earlier, there are six equations with 
seven unknowns related to these problems. Therefore, one of the 
unknowns has to be an independent variable. Consequently, the 
appropriate choice is either (y£) or (<fr) since these two variables are 
related by the equation of conservation of power. Once one of these two 
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variables is fixed, the other variable is then found. The conservation 
of power principle implies that 


(V 


F(4>)d4> = 


C l(y)dy 


or 


F(4>)d* = C I(y)dy 


(4.1) 


(4.2) 


for a differential ray tube in which y Q is the location of central ray 
at the aperture plane. This corresponds to the incident ray coming 
along of the primary feed where ♦ is measured from the positive x- 

axis. From Equation (4.2), one can determine the point by point 
relation between the aperture point y 2 and feed angle $ and thus, the 
rest of the equations for the reflector surfaces are solved by the 
method discussed in Chapter III and Appendix B. 

4.1 Shaping of a Center-Fed Gregorian Reflector Antenna for a Uniform 

Aperture Distribution 

In this section, a center-fed two-dimensional Gregorian reflector 
antenna is used for a starting point of shaping. The geometry of this 
reflector is shown in Figure 4.1 with the following parameters: 

F = 0.6667' 
m 

F = 0.1667' 
c 
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L y = -0.1' 

D = 2.0' 
m 

4» r = 31.42°. 

These parameters are defined in Chapter II. The requirement for shaping 

in this example is to spread the power of the primary feed between -f 

and +$ r into a uniform aperture distribution between y— 1.0' and y=+1.0' 

with the central ray located at y=0.0'. A magnetic line source with 

20 

radiation intensity cos ♦, as shown in Figure 4.2, is used as the 
primary feed. The scattered field patterns are calculated for both the 
unshaped and shaped reflector for the purpose of checking the resulting 
scattered fields after the reflector surfaces are shaped. 

The far-zone scattered fields from the subreflector of the original 
Gregorian reflector antenna calculated at 20 GHz are shown in Figure 4.3 
and the near-zone scattered fields from the parabolic main reflector 
which were calculated at x=2.0' and 20 GHz are shown in Figure 4.4. In 
Figure 4.3, 0 is the angle measured from the negative x-axis and the 
phase center for the subreflector scattered fields is located at focal 
point #2 of the reflector. The dashed line in these figures shows the 
geometrical optics reflected field while the solid line shows the 
scattered fields calculated by the method of physical optics. The 
ripple in the physical optics results is caused by the interaction of 
the geometrical optics field with the edge diffracted field from the 
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Figure 4.3. Far-zone scattered field of the subreflector for the 
conventional Gregorian reflector at 20 GHz. 






sharp terminations of the reflector surface. Consequently, the 
geometrical optics field passes through the average of the ripple. It 
can be seen from the optics field that the main reflector has an 
aperture edge taper of about -14dB. It is thus expected that the 
reflector surfaces have to be shaped significantly in order to achieve a 
uniform amplitude distribution. Also note that the G.O. phase pattern 
of the far-zone reflected fields from the subreflector is constant since 
the reflected field appears to emanate from focal point #2 which is the 
phase center. The phase pattern of the main reflector indicates that 
the aperture field is a plane wave. 

It was mentioned earlier that in order to solve for the shaped 
reflector surfaces, one has to provide a set of appropriate initial 
points. However, if the initial points were not provided correctly, the 
resulting surfaces tend to adjust themselves in order to get the correct 
surfaces. Both of these cases are given next. Example G1 uses a set of 
initial points which do not satisfy the ray condition derived in 
Appendix C; whereas, Example G2 uses appropriate initial points which 
satisfy the ray condition. 

4.1.1 Example G1 

The initial points for shaping are chosen as the intersection 
points of the central ray with the initial subreflector and main 
reflector surfaces, and the initial point on caustic is the second focus 
of the initial elliptic subreflector. The resulting reflector surfaces 
are shown in Figure 4.5. The subreflector surface is replotted in 
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Figure 4.6. The solid line shows the shaped reflectors whereas the 
dashed line shows the initial Gregorian reflector. The initial points 
for the shaping are indicated by small circles in these figures. Also 
shown in Figure 4.5 is the caustic curve for the reflected rays from the 
shaped subreflector. The upper half of the caustic curve corresponds to 
the lower half of the subreflector and the upper half of the main 
reflector; whereas, the lower half of the caustic curve corresponds to 
the remaining reflector surfaces. As is evident in Figure 4.6, the 
choice of the initial points are not quite correct. This happened 
because the power of the primary feed needs to be redistributed, and 
thus, the initially chosen points will not give the correct field at the 
aperture. However, the results of the analysis show where the 
appropriate initial points should be located since the solution tries to 
adjust the incorrect initial points to fit the right surfaces. Thus, 
one can choose the initial points indicated in these results and repeat 
the shaping process. The resulting surfaces which are designated as 
Example G1 are shown in Figures 4.7 and 4.8, and turn out to be 
identical to the previous results except for the initial points. This 
illustrates one of the characteristics of this method of synthesis; 
namely, that incorrect initial points on the reflector surfaces and 
caustic curve will be automatically corrected. 

It can be seen from Figure 4.7 that the caustic of the reflected 
ray from the subreflector does not stay at a fixed point, instead, it 
moves gradually near the beginning of shaping but changes rather rapidly 


64 


Figure 


o 

l/> 



X (FEET) 


4.5. Shaped dual-reflector with inappropriate initial points of 
Example Gl. 
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Figure 4.6. Subreflector and caustics for the shaped dual-reflector 
with inappropriate initial points of Example Gl. 
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Figure 4.7. Improved shaped dual-reflector of Example Gl. 
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Figure 4.8. Subreflector and caustics of the improved shaped dual 


reflector of Example G1 


toward the end of the shaping. The explanation for this behavior is that 
during the shaping process the incident field which emanates from the 
primary source tapers off, and thus, the reflected field from the 
subreflector must also decrease unless the spread factor of the 
reflected field increases. It is the magnitude of this subreflector 
reflected field which determines the magnitude of the aperture field 
since its magnitude does not change after reflection from the main 
reflector. The spread factor for this case is given by 



where is positive and is the distance between the caustic and the 
reflection point on the subreflector; and s^ is the distance between the 
reflection points on the subreflector and main reflector. Thus, in 
order to increase the spread factor the magnitude of has to be 
increased. Consequently, the caustic curve moves away from the 
subreflector surface. Figure 4.9 shows two rays for this reflector. 

The scattered fields for Example G1 which are calculated at 20 GHz 
are given in Figures 4,10 and 4.11 for the far-zone subreflector and 
near-zone main reflector patterns, respectively. Both of these patterns 
are different from the conventional Gregorian reflector as they should 
be. The resulting near-zone scattered fields from the main reflector 
indicate that a uniform plane wave with a uniform amplitude distribution 
is nearly obtained. Also notice that the phase pattern of the 
subreflector geometrical optics scattered fields is not a constant as in 
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9. Ray geometries for the shaped reflectors of Example Gl. 
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the conventional Gregorian reflector case since the caustic of the 
shaped subreflector reflected field is not a fixed point. 

4.1.2 Example G2 

The initial points on the subreflector and main reflector are 
determined by the intersection of the central ray with the original 
subreflector and main reflector of the Gregorian antenna in this 
example. However, instead of using the second focus of the elliptic 
subreflector as the initial point for the caustic curve, the ray 
condition discussed in Appendix C is used to determine the initial point 
for the caustic. The resulting surfaces which are designated as Example 
G2 are given in Figures 4.12 and 4.13 for the overall system and the 
subreflector, respectively. Figure 4.14 shows two rays for the shaped 
reflectors. The far-zone scattered fields of the shaped subreflector 
and the near-zone scattered fields of the shaped main reflector 
calculated at 20 GHz are shown in Figures 4.15 and 4.16. Comparing 
these results with the ones of example Gl, it is found that the 
subreflector scattered patterns are not quite the same, but the 
scattered patterns of the main reflector are identical in magnitude. 

This shows the non-uniqueness of the shaping solution when the initial 
points for shaping are different. The major difference in the reflector 
surfaces between Examples Gl and G2 is that the subreflector surface and 
the caustic curve are shifted although the overall surfaces are only 
slightly different. 
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Figure 4.12. Shaped dual-reflector of Examp 
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Subreflector and caustics for the shaped dual-reflector of 


Example G2 
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4.2. Shaping of an Offset Reflector Antenna (Example G3) 

In the previous section, a center-fed symmetric Gregorian antenna 

was shaped to generate a uniform plane wave with a uniform aperture 

distribution in which the aperture was designated between y— 1' to 

y-+l'. In this section, an offset fed Gregorian reflector is used as 

another example to illustrate that the same shaping process can be 

applied to an offset dual-reflector system. This example is designated 

as Example G3. The geometry of the initial offset Gregorian reflector 

is as follows: F =9.5', F =8.5', L =-1.5' and is shown in Figure 4.17. 
m c v 

The central ray of the main reflector is designated at y=8.5'. The 
initial points on the subreflector and main reflector surfaces are 
chosen as the intersection points between this central ray and the 
initial Gregorian antenna. Consequently, the primary feed axis must be 
tilted by an angle of -7.68°. The primary source is a magnetic line 
source with a radiation intensity cos^% with respect to the feed axis. 
The feed pattern is shown in Figure 4.18. A plane wave with a uniform 
amplitude distribution is to be obtained between y=-1.0' and y*18.0' by 
shaping the given offset Gregorian reflector antenna. This results in 
most of the power of the primary source being contained between 4*-8.5° 
and 4=8.5° with respect to the tilted feed axis which illuminates the 
subreflector. The feed pattern illuminates the edge of the subreflector 
at about 9.6 dB below the peak of the feed pattern. 

The far-zone scattered field from the original subreflector 
calculated at 3 GHz is shown in Figure 4.19 for O°<0<18O° where 0 is 
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Figure 4.19. Far-zone scattered fields from subreflector of the offset 
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Gregorian reflector with a cos zut 9 magnetic line source 
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measured from the negative x-axis. The near-zone scattered field of the 
original main reflector calculated in the x=20.0' plane and at 3 GHz is 
shown in Figure 4.20. It can be seen from Figure 4.20 that the original 
offset Gregorian reflector has an aperture edge taper of about 10 dB. 

The initial point on the caustic curve is determined by the ray 
condition, and the resulting shaped reflectors are given in Figures 4.21 
and 4.22. Figure 4.23 shows two rays for the shaped reflectors. The 
far-zone scattered fields from the shaped subreflector and the near-zone 
scattered fields from the shaped main reflector calculated at 3 GHz are 
shown in Figures 4.24 and 4.25, respectively. It is obvious that the 
requirement of obtaining a uniform aperture distribution is 
approximately achieved except for the ripple which is caused by the 
interaction of the edge diffracted field with the reflected field. 

4.3. Shaping of a Center-Fed Cassegrain Reflector Antenna for Uniform 
Aperture Distribution (Example Cl) 

An example of shaping a center-fed Cassegrain reflector is 
discussed in this section. The geometry of the original antenna is 
given in Figure 4.26 with the following parameters: 

F =0.6667' 
m 

F =0.5667' 
c 

L =0.1' 
v 
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Figure 4.23. Ray geometries for the shaped reflectors of Example G3. 
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Figure 4.24. Far-zone scattered fields from the subreflector of the 

shaped offset dual-reflector with a cos^% magnetic line 
source illumination at 3 GHz. 
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Figure 4.26. A conventional Cassegrain antenna. 
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20 

A cos^ ♦ power pattern magnetic line source is used as the primary 
source. The incident field at the edge of the subreflector is 4.5 dB 
down from the maximum. A uniform plane wave with a uniform amplitude is 
to be generated between y=-l' and y=l' in which the central ray is 
defined at y*0. The initial points for the shaped reflectors are chosen 
at the intersection points between the central ray and the original 
reflectors. The initial point on the caustic curve is thus determined. 
The resulting shaped reflectors are shown in Figures 4.27 and 4.28 with 
the original reflectors indicated by dashed lines and the shaped 
reflectors by solid ones. The upper half of the caustic curve 
corresponds to the lower half of the subreflector and main reflector; 
whereas, the lower half of the caustic curve corresponds to the 
remaining surfaces. The behaviour of the caustic curve can be 
explained in the same way as for the Gregorian case which is given in 
Example Gl. The caustic curve moves away from the subreflector surface 
in order to increase the reflected caustic distance p c , which is given 
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Figure 4.27. Shaped dual-reflector with convex subreflector of Example 
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Figure 4.28. Subreflector with the convex surface of shaped dual 


reflector Example Cl 
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Note that p increases in order to increase the reflected field and thus 
c 

compensate for the taper of the primary source pattern. 

The scattered fields of the original subreflector and main 
reflector calculated at 20 GHz are shown in Figures 4.29 and 4.30. The 
scattered fields of the shaped reflectors are shown in Figures 4.31 and 

4.32. Two ray geometries for the shaped reflectors are shown in Figure 

4.33. 


4.4. Shaping of a Center-Fed Gregorian Reflector Antenna for Non- 

Uniform Aperture Amplitude Distribution (Example G4) 

In the previous examples, the purpose of shaping was to generate a 

uniform amplitude distribution across a given aperture. In this 

section, an example is used to illustrate that the same shaping process 

can be used to generate a non-uniform amplitude distribution. The 

Gregorian reflector antenna of example G1 and a magnetic line source 
20 

with cos $ power pattern are used to generate the following aperture 
distribution: 


( 1 


i(y) 




|y I <o.5' 

l'>|y|>0.5' 


This aperture distribution which is shown in Figure 4.34 has a -lOdB 
edge taper at the edge of the main reflector; i.e., y=±1.0'. The 
central ray is thus at y=0' . 


96 


The resulting shaped reflectors are given in Figures 4.35 and 4.36. 
It can be seen from these results that the caustic curve has a 
discontinuity. The explanation for this behavior is that since the 
aperture field is a function of 



vhere is the distance between caustic and the reflection point on the 
subreflector, the shaping process tried to make smaller in order to 
generate a rapidly decreasing subreflector reflected field and aperture 
distribution. Figure 4.37 shows two rays for the shaped reflectors. 

The far-zone scattered fields of the shaped subreflector are given 
in Figure 4.38; while, the near-zone scattered fields calculated in the 
x-2.0' plane of the main reflector are given in Figure 4.39. These 
scattered fields are calculated at 20 GHz. One can see that the 
amplitude of the ripple on the near-zone scattered fields is smaller 
than for the uniform aperture distribution case (i.e., Examples G1 and 
G2) since the incident field at the main reflector edges is smaller, the 
edge diffraction is not as significant. 

4.5 Summary 

Several examples are given in this chapter to demonstrate the 
shaping process discussed in Chapter III. Scattered field results are 
calculated for each case. These results indicate that the shaped 
reflector can provide considerable improvement in obtaining the desired 
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Figure 4.34. Desired near-zone aperture distribution of Example G4. 
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Figure 4.36. Shaped subreflector and caustic of Example G4 
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Figure 4.38. Far-zone scattered fields from shaped subreflector of 
Example G4 at 20 GHz. 
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aperture distributions, as compared to the conventional Gregorian and 
Cassegrain reflectors. The geometrical nature of the synthesis method 
provides more insight into the way that the shaped reflectors work than 
do the methods which depend on mathematical solutions to a coupled set 
of differential equations. A major advantage is that one can determine 
much about the practicality of the resulting reflector surfaces by 
examining the caustic curve of the shaped subreflector. A good example 
is the discontinuity in the caustic curve shown in Figure 4.36 for 
Example G4. The kinks of the caustic curve indicate that the slope of 
the surface curvature of the shaped subreflector is discontinuous which 
cannot easily be seen by looking directly at the subreflector. 

It should also be noted that the amplitude of the resulting 
secondary aperture distribution is controlled by the amplitude of the 
subreflector reflected field at the main reflector surface. The main 
reflector is used primarily to adjust the phase of the reflected field 
such that the resulting aperture distribution is uniform in phase. 

All the pattern results shown in this chapter include the 
diffracted fields from the edge of the reflectors. For compact range 
applications, this will not provide satisfactory performance. In the 
next chapter, a physical optics analysis is used to determine the 
performance of the shaped main reflector with a blended rolled edge 
attachment in order to reduce the edge diffracted fields. 
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CHAPTER V 


RO LLED EDGE MODIFICATION FOR THE SHAPED MAIN REFLECTOR 
5.1 Introduction 

For far-zone scattering and antenna measurements, one is required 
to have a plane wave with a uniform amplitude illuminate either the 
scatterer or antenna under test. This is often achieved by using an 
outdoor far-zone range to simulate the plane wave illumination. There 
are some serious limitations associated with outdoor measurements, such 
as weather, interference, etc., which limit the accuracy of measurement. 
Consequently, indoor measurement facilities are an alternative to 
outdoor measurements although the uniform plane wave illumination is 
still a requirement. The creation of a uniform plane wave by a focussed 
parabolic reflector has been an attractive approach for indoor 
measurement facilities. Conventional parabolic reflector antennas have 
been used in so-called compact range measurement facilities [22]. 

Compact ranges are indoor facilities which provide the capability to 
measure antenna pattern and scattering by simulating a plane wave 
illumination on the test antenna or scatterer. There are basically two 
drawbacks associated with the use of reflectors in compact range 
applications. The first is that the plane wave reflected from the 
parabola is not uniform in amplitude, as shown in Figure 5.1, due to the 
illumination taper of the feed-horn and space-attenuation effects. 
Consequently, the size of the available measurement volume is limited, 
thus limiting the size of the test antenna or scatterer. The second 
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Figure 5.1. Aperture edge taper due to feed-horn pattern and space 
attenuation. 
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drawback is that the abrupt termination of the reflector surface which 
creates a very strong edge diffracted field as shown in Figure 5.2. 

This diffracted field interferes with the plane wave and results in 
magnitude and phase variations of the wave illuminating the test antenna 
or scatterer. The measured results thus may not be accurate. 

As pointed out in [2] the nonuniformity of the near-zone 
distribution can be improved by using an offset reflector system such as 
a Gregorian antenna with the primary feed tilted properly. It is also 
conceivable that by properly shaping a dual-reflector antenna, the edge 

taper problem can also be improved. 

To solve the problem of edge diffractions caused by the sharp 
termination of the reflector surface, one can make the reflector size as 
large as possible so that the edge diffracted fields are not very strong 
in the region where the antenna or scatterer under test is located so 
that the resulting variations in the plane wave are not as significant. 
However, this approach is very costly because it is not practical to 
build such a huge reflector. The other way of reducing the edge 
diffracted fields is to modify the sharp termination of the reflector by 
adding rolled edges [23], as shown in Figure 5.3. The idea of a rolled 
edge modification has been extensively studied by Burnside et al. 

[1,2, 3, 4]. The addition of rolled edges to the reflector surface 
introduces reflected fields from both the rolled edges and the 
reflector. These two reflected fields are not continuous and there are 
still diffracted fields which emanate from the junction of the rolled 
edge and the reflector. However, this discontinuity in the reflected 
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Figure 5.2. Variations of aperture field due to the interference 
between reflected and diffracted fields. 



JUNCTION 




fields is much smaller than for the reflector with a sharp termination 
since in the latter case, the GO reflected field vanishes outside the 
shadow boundary of the reflector. Consequently, the variations in the 
aperture fields are not as significant when the reflector is modified 
with a rolled edge. 

In this chapter, the rolled edge idea is implemented in the 
analyses of the shaped main reflector in order to show the improvement 
in the performance of the shaped dual-reflector for compact range 
applications. Some examples are presented to validate this concept. In 
calculating the scattered fields from the shaped main reflector and the 
rolled edges, the surface currents are assumed to exist only in the 
region which is illuminated by the reflected fields from the 
subreflector and vanish elsewhere, as shown in Figure 5. A. This 
assumption creates false diffracted fields from the points where the PO 
currents vanish. The false diffracted fields also interfere with the 
true scattered fields from the reflector surface. The method for 
correcting this problem has been studied by Gupta and Burnside [3] and 
is also discussed and used in this chapter. 


5.2 Rolled Edge Attachment to Main Reflector 

As was discussed earlier, the purpose of adding a rolled edge to 
the main reflector termination is to create a smooth transition of the 
geomtrical optics reflected fields from the reflector to the rolled edge 
in such a way that the diffraction from the junction of the reflector 
and the rolled edge is insignificant. Using this approach, the 
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Figure 5.4. Induced currents on main reflector with rolled edge 


modification. 



resulting incident wave in the measurement zone has very small ripple. 
One other requirement of the rolled edge is that the reflected fields 
from the rolled edge do not enter the measurement volume to interfere 
with the reflected fields from the reflector. The detailed analysis of 
a blended rolled edge attachment to the parabolic reflector can be found 
in [2]. In this section, the rolled edge is used to modify the edge of 
the shaped main reflector. The basic concept is still the same as for 
the parabolic reflector case. 

5.2.1 Elliptic Rolled Edge 

Consider an elliptic rolled edge attachment to the shaped main 
reflector as shown in Figure 5.5. The shaped main reflector is 
terminated at point P with surface normal n and the ellipse is attached 
at F in such a way that the tangents of the main reflector and the 
ellipse are continuous. The two semi-axes of the ellipse are given by 
a @ and b @ . Consider a local coordinate system for the ellipse in which 
the coordinate origin is at P, and x e , y g are the two axes such that 
x *n, y -t. The parametric equation for the ellipse is given by 

x g = a e (cosv-l) (5.1) 

and 

y = b sinv (5.2) 

e e 

for 0<v<v 
- - m 
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Note that v is the angle measured from the x g axis, and \> m is the angle 
where the ellipse stops. In terms of the original x-y coordinate system 
the surface of the ellipse can be expressed by 


c ... = a (cos\>-l)cos0.-b sinvsinG. +x. 

ellipse e v t e t j 


and 


( 5 . 3 ) 


y ... = a (cosv-l)sin©. + b sinvcos©. +y . 

■'ellipse e v t e t j 


( 5 . 4 ) 


in which 0 t is the angle between n and the positive x-axis. This 
defines the elliptic rolled edge attachment to one edge of the main 
reflector. A similar definition can be used for the other edge. Since 
the tangent at the junction of main reflector and ellipse is continuous, 
the diffracted fields from this junction will be much smaller than for a 
sharp termination of the main reflector. However, as seen from Figure 
5.5, the main reflector is a concave surface; whereas, the ellipse is a 
convex surface. The radius of curvature at P changes sign between the 
main reflector and the ellipse and thus is obviously discontinuous. 

This results in a discontinuity of the G.O. reflected field between the 
main reflector and ellipse. Consequently, diffracted fields emanate 
from this junction in order to compensate for the discontinuity in the 
reflected fields. In order to further reduce this diffracted field, a 
blended rolled edge which is created by blending the main reflector with 
the ellipse by some functions is used to improve the performance of the 
rolled edge. This is discussed in the next section. 
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5.2.2. Blended Rolled Edge 

As discussed In the previous section the change of surface from the 
main reflector to the ellipse creates diffracted fields from the 
junction. A gradual change of the surface from the main reflector to 
the ellipse can further reduce the diffracted fields from the junction. 
This is done by the so-called "blended" rolled edge modification. The 
elliptic rolled edge discussed in the previous section is transformed 
into a blended rolled edge by the following relationship: 


W’> ■ 'empse< v)b<v > + W v)<1 - b( ' ,)) <5 ' 5) 

in which / ellipse ( v > is the equation of the ellipse, / main (^) is the 
equation of the main reflector extended beyond the junction and b(v) is 
the blending function such that 0 < b(v) <1 and b(0) - 0, b(v m ) - 1. 
The idea of blending is shown in Figure 5.6. Note the the blended edge 
is attached in the same way as the elliptic edge so that the surface at 
P is smooth and continuous. 

The radius of curvature of the surface at the junction is also 
continuous since the blended edge at the junction is actually the main 

reflector because b(0)=0, and *blend <0 > ■ W°>' There are many 
possible blending functions b(v) such as 


1) b(v) = 

m 


linear blending 
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2) b(v) 


: square blending 



3) b(v) ■ i (l - cos — 1 : cosine blending, and 

L 1 v m J 

- cos : cosine square blending 

All of these functions satisfy the condition that 0<b(v)<l, b(0)=0 
and b(v m )=l. A detailed analysis of these blending functions can be 
found in [2]. 

In order to generate the blended edge, it is necessary to specify 
the portions of the pure ellipse and the extended main reflector surface 
to be used. This is done by appropriately choosing the values of the 
parameters a , b , v , and s as shown in Figures 5.6 and 5.7. The 
portion of the ellipse is defined by the maximum elliptic angle v . The 
portion of the main reflector is defined by s m , the maximum distance 
along the tangent of the surface at the junction point P. The point on 
the ellipse which corresponds to a given point (x 2 ,y 2 ) on the extended 
main reflector is defined as follows. The parameter s along the tangent 
t is calculated first by 


A) b(v) * \ 


s 




where m^ is the slope of t. 


(5.6) 
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Figure 5.7. Definition of parameters for blending. 
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Then, the elliptic angle v which defines a point on the ellipse is 
related to a point on the extended reflector surface through the 
parameter s by 


v = — v 
s m 
m 


(5.7) 


From v, the point on the ellipse corresponding to x 2 ’ y 2 can be 
calculated by Equations (5.3) and (5.4). Finally, the blended rolled 
edge is obtained from Equation (5.5) and is given by 


x b * (l-b(v)]x 2 + b(v)ta e (cosv-l)cos0 t - b e sinvsine t +Xj ] (5.8) 


y b = [l-b(v)]y 2 + b(v)[a e (cosv-l)sine t + b e sinvcos6 t +y ; . ] (5.9) 


Although the slope and surface radius of curvature at the junction 
between the main reflector and blended rolled edge are continuous, there 
are higher order derivatives of the surface which are not continuous at 
the junction and thus, some diffractions emanate from the junction. The 
diffraction coefficient for such diffraction is not available so that 
only the method of Physical Optics can be used to analyze the 
performance of the reflector surface with a blended rolled edge. 

However, as mentioned earlier, the induced surface currents in the 
physical optics integral are assumed zero on the shadow side of the 
surface although this is not true in reality. A method for correcting 
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this problem has been investigated in [3] and is discussed in the 
following section. 

5.3 End-Point Correction of Physical Optics 

The scattered fields calculated by the method of physical optics 
for a two-dimensional scatterer are given by Equations (2.21) and 
(2.22). For a magnetic line source illuminating on the scatterer, the 
scattered magnetic fields are given by 



in which ? is the induced physical optics surface current on the 
* 

scatterer, p is the unit vector in the observation direction from the 
current and p is the distance between the observation point and the 
current. The physical optics current, ^ , exists over the region 

directly illuminated by the line source and is assumed to vanish on the 
shadow side of the scatterer. The integration is thus performed only 
over the lit side of the scatterer. By applying Equation (5.10) to the 
shaped main reflector with a blended rolled edge termination on both 
ends of the main reflector section, one can compute the ordinary P0 
scattered fields from the entire reflector. The scattered fields 
calculated this way include the following contributions: 
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1. The stationary phase or reflection point contribution 
which is equivalent to the geometrical optics reflected 
fields, 

2. The two end-point contributions which result from the 
sharp termination of the surface currents at each end 
shadow boundary, and 

3. The diffracted fields from the two junctions between the 
main reflector and the blended rolled edge. 

These contributions are illustrated in Figure 5.8. The end-point 
contributions are incorrect contributions because surface currents still 
exist on the shadow side of the rolled edge. It is stated in [3] that 
these end-point contributions do not give the correct creeping wave of 
the curved surface. Consequently, the end-point contributions have to 
be removed from the physical optics integral. This can be accomplished 
by evaluating the physical optics integral asymptotically. 

Consider the integral given by 


I 


r p(i)* w(i) di . 


(5.11) 


The asymptotic evaluation of this integral [3] results in 


I ~ 


2n 


T 


, k |f "(l g ) | 


jk/(l ) jn/4 sgn[r<l )] 

F(l g )e S e 


00 



n=0 


<-D n F n 

jk/' 


b 


e jkf(l) 


a 


(5.12) 
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Figure 5.8. Illustration of three contributions from the physical 
optics integration. 
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where 1 is the stationary phase point and is assumed away from the two 
s 

end-points l=a,b. Note that F n is defined as 

F „ - 3i fem] <5 - 13) 

with 


F q ( 1) = F(l) . (5.14) 

The first term in Equation (5.12) is the specular reflected field and 
the second term represents the two end-point contributions. Since 
Equation (5.12) is an asymptotic representation of (5.11), there are 
some higher order terms which have been neglected. One of these terms 
is the diffraction from the junction of the main reflector and the 
blended rolled edge. Consequently, one obtains the desired result by 
subtracting these end-point contributions from the PO result such that 
the scattered field is given by 


H 


s 


’W k « 1 >41 


“ (-l) n F 


E 

n=0 


jk/' 


_n_ fi jk/(l) 
( 1 ) 


(5.15) 


Thus, by numerically integrating Equation (5.11) and subtracting the end 
point contributions, one obtains the geometrical optics reflected field 
plus the diffracted fields from the junction of the main reflector and 
rolled edge. For surfaces which can be described analytically such as a 
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parabola, the end-point contributions can be calculated analytically 
[2]. In this study, the surface is generated numerically, and thus, the 
end-point contributions are computed numerically. The rest of this 
section discusses the computation of the physical optics currents on the 
shaped main reflector by geometrical optics for the cases where the 
shaped subreflector is either concave or convex. 

In formulating the surface equations for the shaped dual-reflector, 
it is assumed that the surface points (x^,y^) and (X 2 »y 2 ) are located on 
either a Cassegrain or Gregorian reflector with the caustic of the 
subreflector reflected field at (x c ,y c ). Thus, the equations for 
calculating the geometrical optics reflected field of a Cassegrain or 
Gregorian reflector are used to calculate the reflected field of the 
shaped reflector. For the case of a shaped dual-reflector with a 
concave subreflector, as shown in Figure 5.9, the surface points (x^,y^) 
and (X 2 »y 2 ) are located on a particular Gregorian antenna with the 
caustic of the subreflector reflected field located at (x c ,y c ). From 
Equation (C.9), the reflected field at (* 2 ^ 2 ^ from the subreflector can 
be calculated by 


H J< x 2 »y 2 > - H z (s i> = H z (x i»yi> 


-p„ 


-jks 1 


- p c +s l 


(5.16) 


in which H^(x^,y^) is the incident magnetic field at (x^,y^) and is 
the distance between (x^,y^) and ( x c »y c )> and s i * s distance between 
(Xj^f) and (x 2 ,y 2 )* The incident magnetic field H^(x^,y^) is given by 
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(5.17) 


H z< x i»yi> = J p ^i > 1 


-jkpj 

e 



where P(^) is the far-zone power pattern of the magnetic line source at 
angle The phase of the far-zone field pattern of the line source is 

assumed constant and suppressed. Thus, the induced physical optics 
surface current on the main reflector is given by 

^po = 2 n x ^(x 2 ,y 2 ) = 2 n x i r (s 1 ) = -12 H*(x 2 ,y 2 ) 
or 

3 = i J, (5.18) 

po 1 


where 


J 


1 


-2|PH^T 


-jk Pi 

e 



-P c -jks 1 

-p + s. e 
i c 1 


(5.19) 


Consequently, the PO scattered fields from the main reflector are 
calculated by 


H“(x,y) = 


l4n 


J n/4 


m * 


-jkp 


z • (1 x p)J x 
O ip 


dl 


or 
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H®(x,y) 



-jk p i 

• (1 x p)(-2) 

i*l 


p c "j ks i e -j k P 


^-p c+ s i 




dl . 


(5.20) 


For the surface current on the blended edge point (x^y^), one has to 
numerically find the corresponding incident ray from the subreflector; 
i.e., corresponding to (x^y^ and (x c ,y c ). By comparing Equations 
(5.20) and (5.11), it is found that 


F(l) = 2 


^2 e JnM G • (i x p)) -IPC+j)’ 




-p„ 


■ p c +s l 


(5.21) 


and 


f (1) = - (pj+Sj^+p). 


(5.22) 


Similarly, for the case of the convex subreflector, the incident 
field at the main reflector is given by 


H z (x 2^2> = H z (s l> = H z (x l’ y l> 


-jks. 


s l +p c 


(5.23) 


where (x^y^ and (x 2 ,y 2 ) are the points on the particular Cassegrain 
reflector with subreflector caustic located at ( x c »y c > as shown in 
Figure 5.10. The surface currents are then given by 
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with 


J x = -2 4P(^) 


i e 


-jkp. 

If 


-jks. 


s l +p c 


(5.24) 


Finally, the PO scattered field from the main reflector is calculated by 


H®(x,y) 


Jn/4 
4 ii 6 


Urn 


z*(i x p) (— 2 ) Jp (^T 


-jkpi 

e 



-jks 1 


. P +S., 

■i y c 1 


-jkp 

Hr 


dl 


(5.25) 


This gives 


F(l) = 2 


\ jn/4 
4 n 6 


(i x P ) 4 p<*JF 

|pTIp 


p c +s l 


(5.26) 


with f ( 1) given in Equation (5.22). 

Examples are given in the next section to illustrate the 
combination of a shaped dual-reflector with blended rolled edge 
terminations to the main reflector. The scattered fields are calculated 
for the main reflector with a blended edge. The results of the PO 
integration, the PO integration with the first order end-point 
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Figure 5.10. Surface current calculation for the case of shaped dual 
reflector with convex subreflector. 
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corrections, and the PO integration with the first and second orders 
end-point corrections are, shown to illustrate the effects of the false 
end-point contributions of the ordinary physical optics integral. 
Normally, the rest of the higher order end-point contributions are very 
small and thus are neglected in the following examples. 

5.4. Examples of Shaped Dual-Reflectors with Blended Rolled Edge 
Terminations on the Main Reflector 
5.4.1. Example BL1 

In this example, the Gregorian reflector antenna given in Examples 

G1 and G2 of Chapter IV is used as the starting point for the shaping, 

and the cosine-squared blending function is used to generate the blended 

rolled edges which are attached to the main reflector at y=0.5' and 

y=-0.5' . The parameters for the blended edges are given by a^O.2', 

b *0.5' , s =1.2' and v = 90°. Since the blended rolled edges require 
e m m 

an extension of the shaped main reflector, the shaping of the main 

reflector is extended to y=1.3' and y=-1.3'. The aperture field within 

this region is assumed to be uniform. A magnetic line source with a 

20 

power pattern of cos <f> is used as the primary source. The tilt angle 
of the feed axis is 0° and the power between 4=-40° and +=40° is 
redistributed uniformly over the region between y=-1.3' and y=1.3'. 

The resulting shaped reflectors are shown in Figure 5.11. Note that 
these surfaces are different from those of Example G2 which are shown in 
Figure 4.12. This is due to the difference in the designated aperture 
size and the amount of the feed power which is redistributed. As seen 
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from Figure 5.11, the caustic curve intersects the subreflector surface 
so that one might think that some of the reflected field from the 
subreflector will be blocked by the subreflector itself. However, it is 
found by tracing the rays that only the reflected fields from the 
central part of the subreflector illuminate the blended main reflector. 
This is shown in Figure 5.12. Therefore, the scattered fields from the 
blended main reflector can still be calculated to validate the 
performance of the blended rolled edges. The near-zone scattered fields 
calculated by the method of PO at 20 GHz are shown in Figure 5.13. The 
scattered fields with the first order end-point corrections are shown in 
Figure 5.14, and the scattered fields with both first and second order 
end-point corrections are shown in Figure 5.15. One can see the effects 
of the false end-point contributions by comparing these results. By 
taking out these false contributions, the scattered fields within the 
unblended section of the whole main reflector, i.e., y=o' to y=0.5' in 
Figure 5.15, are very smooth and uniform. This is a significant 
improvement in the ripple of the scattered fields for compact range 
applications compared to the one without blended edges. The slight 
ripple shown in Figure 5.15 is caused mainly by the diffraction from the 
junction of the unblended and blended sections of the main reflector. 
Also note that if the "true" scattered fields from the subreflector 
(i.e., not just the reflected fields from the subreflector) are used to 
calculate the surface currents on the main reflector, there will be 
slowly varying ripple in addition to the scattered fields from the main 
reflector. 
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Figure 5.12. Section of subreflector from which the reflected fields 


illuminate the blended main reflector in example BL-1. 



5.4.2. Exaaple BL2 

The same Gregorian reflector antenna used in the previous example 
is also used here as the starting point for shaping. However, instead 
of a uniform aperture distribution, the distribution given in Figure 
4.34 of Chapter IV is used and the blended rolled edges are attached to 
the main reflector at 7=10.5'. The parameters associated with blended 
rolled edges, feed pattern, and feed power to be redistributed over the 
designated aperture, are all the same as in the previous example. 

The resulting shaped reflectors with blended rolled edge 
terminations to the main reflector are shown in Figure 5.16. Again, the 
section of the subreflector from which the reflected fields illuminate 
the blended main reflector is replotted in Figure 5.17. 

The near-zone scattered fields calculated at 20 GHz from the 
blended main reflector are given in Figures 5.18 through 5.20 for the 
PO, the PO with first order end-point corrections, and the P0 with first 
and second order end-point corrections, respectively. A very smooth 
result is obtained between y*0' and y=0.5' when the first and second 
order end-point corrections are used. However, if one compares the 
pattern shown in Figures 5.15 and 5.20 where the first one is obtained 
by shaping reflectors to generate a uniform aperture distribution while 
the second one is obtained by shaping the reflectors to generate a non- 
uniform aperture distribution, the ripple in Figure 5.20 between y=0.5' 
and y*0.7' is higher than for the one in Figure 5.15. This shows that a 
potential problem can occur if one tries to shape the reflectors to 
generate an aperture distribution which is constant over one region and 
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of Example BL-1 at 20 GHz using P0 alone. 
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Figure 5.14. Near-zone scattered fields from the blended main reflector 
of example BL-1 using PO with first order end-point 
corrections at 20 GHz. 
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Figure 5.15. Near-zone scattered fields from blended main reflector of 
example BL-1 using PO with first and second order end- 
point corrections at 20 GHz. 
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Figure 5.16. Shaped dual reflectors with blended rolled edge 
terminations for example BL-2. 
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Figure 5.17. Section of subreflector from which the reflected fields 


illuminate the blended main reflector in example BL-2 
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Figure 5.18. Near-zone scattered fields from blended main reflector of 


example BL-2 using PO at 20 GHz. 
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of example BL-2 using PO with first order end-point 
corrections at 20 GHz. 
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of example BL-2 using PO with first and second order end- 
point corrections at 20 GHz. 
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drops off rapidly over the other region. This can cause some 
diffractions which will create significant ripple in the desired 
aperture distribution. 

5.4.3. Example BL3 

The offset Gregorian reflector used in Example G3 of Chapter IV is 
used as the starting point for shaping in this example. A uniform 
aperture distribution over the region between y=-2.6' and y=18.8' is 
used to generate the shaped dual reflector. The blended rolled edge 
begins at y=5.5' for the lower edge and at y=11.5' for the upper one. 

A magnetic line source with a power pattern of cos t is used as the 
primary source. The tilt angle of the feed axis is -7.7° so that the 
central ray is located at y=8.5'. The power between -9° and 10° with 
respect to the tilted feed axis is redistributed uniformly over y=-2.6' 
to y*18.8'. The parameters for the blended rolled edges are as follows: 

a =0.78', b =4.5', s =8.5', and v = 90° 
e e m in 

and the cosine-squared blending function is used. The resulting 
reflector surfaces are shown in Figure 5.21. 

The near-zone scattered fields calculated at 3 GHz for this new 
reflector system are shown in Figures 5.22 through 5.24 for the P0, the 
P0 with first order end-point corrections, and the P0 with both the 
first and second order end-point corrections, respectively. Similar 
results calculated at 10 GHz are shown in Figures 5.25 through 5.27. 
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Figure 5.22. Near-zone scattered fields from the blended main reflector 
of example BL-3 at 3 GHz using PO. 
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Figure 5.23. Near-zone scattered fields from the blended main reflector 


of example BL-3 at 3 GHz using PO with first order end- 
point corrections. 
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Figure 5.24. Near-zone scattered fields from the blended main reflector 
of example BL-3 at 3 GHz using PO with first and second 
order end-point corrections. 
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Figure 5.25. Near-zone scattered fields from blended main reflector of 
example BL-3 at 10 GHz using P0. 


151 


o 



Figure 5.26. Near-zone scattered fields from blended main reflector of 
example BL-3 at 10 GHz using P0 with first order end-point 
corrections. 
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Figure 5.27. Near-zone scattered fields from the blended main reflector 
of example BL-3 at 10 GHz using P0 with first and second 
order end-point corrections. 
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The very small difference between Figure 5.26 and Figure 5.27 reveals 
that the second order end-point corrections are insignificant in this 
case. The ripple shown in Figure 5.27 results from the numerical 
integration errors associated with the PO solution. 


5.5. Summary 

The addition of rolled edges to the main reflector reduces the edge 
diffracted fields and consequently, reduces the ripple in the near-zone 
scattered fields, as compared to the reflector with a sharp termination. 
This is very helpful in improving the performance of the compact range 
reflector. Several examples are used to verify this concept. The 
physical optics method is used to calculate the scattered field from the 
main reflector with blended rolled edges. Corrections are made for the 
false end-point contributions which result from the physical optics 
approximation. The resulting scattered fields of these examples show 
that very good performance of the reflectors are achieved. 
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CHAPTER VI 


SHAPING OF THREE-DIMENSIONAL CIRCULARLY SYMMETRIC DUAL-REFLECTOR 

ANTENNAS 


Discussions in the previous chapters have concentrated on the 
shaping of two-dimensional dual-reflector antennas. However, as was 
mentioned in Chapter III, the same approach can be applied to three- 
dimensional dual-reflector antennas with circular symmetry and 
circularly symmetric feed pattern illuminations. The approach for 
circularly symmetric shaped reflectors is discussed briefly in this 
chapter. 

6.1. Method of Shaping 

A three-dimensional circularly symmetric surface is a body of 

revolution which is generated by rotating a generating curve about its 

axis of symmetry. For example, a three-dimensional circularly symmetric 

Cassegrain antenna can be generated by rotating the curves shown in 

Figure 2.6 about the x-axis, the axis of symmetry. By using the 

conventional cylindrical coordinates (p,+,z), one can replace the x and 

y axes in Figure 2.6 by z and p axes, respectively. Similarly, if one 

defines the generating curves for a shaped dual-reflector for 0<p<p , 

max 

as shown in Figure 6.1, a three-dimensional circularly symmetric shaped 
dual-reflector can be found by rotating these curves about the z-axis, 
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provided that the given primary source has a circularly symmetric 
pattern, and the desired aperture distribution is also circularly 
symmetric. 


One major difference between the two-dimensional case and the 
three-dimensional one is that the formula for calculating the 
geometrical optics field is different. As discussed in Chapter II the 
geometrical optics field is proportional to a spread factor which 
involves the radius of curvature of the geometrical optics wavefront. 

For the two-dimensional case, the ray tube is two-dimensional, and the 
wavefront has only one radius of curvature. For the three-dimensional 
case, the ray tube is three-dimensional, and the wavefront has two radii 
of curvature. Consequently, one has to integrate the power pattern of 
the geometrical optics field over a surface in order to calculate the 
power contained in a three-dimensional ray tube. 

Consider the shaping of a three-dimensional dual-reflector antenna 
to generate a symmetric near field aperture distribution CI(p) when a 
primary source with a symmetric radiation pattern F(0) is given. 
Conservation of power states that 



F(0) sin9 dGd* = 



CI(p)pdpd* 


( 6 . 1 ) 


in which 6 is the maximum illumination angle at the subreflector 
max 

edge, and p is the maximum radius of the secondary aperture as 
max 

defined in Figure 6.2. From Equation (6.1), the constant C is found as 
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max 


F(0)sin0d0 


( 6 . 2 ) 


max 


I( p) pdp 


Therefore, the relationship between 0 and p is given by 


f0 


F(0) sin0d0 


0 

max 


F(0)sin0d0 


max 


I(p)pdp 


i: 


I(p)pdp 


(6.3) 


It is obvious that the tilt angle of the feed axis is 0° in the 
circularly symmetric case. Equation (6.3) shows that one can use the 
identical shaping process given in Chapter III to create generating 
curves for three-dimensional circularly symmetric dual-reflector 
antennas by using this equation as the appropriate equation for 
conservation for power. The coordinates x and y used for the two- 
dimensional case are changed to z and p for the three-dimensional case. 

The initial points for the shaping process have to be chosen as in 
the two-dimensional case. The condition which these initial points have 
to satisfy is given in Appendix C. Normally, these initial points are 
located on the axis of symmetry. Also the caustic curve in the two- 
dimensional case becomes a caustic surface in the three-dimensional 
case. 


159 



Three examples are given in the next sections to illustrate the 
shaping of three-dimensional circularly symmetric dual-reflector 
antennas. The results of the first two examples are compared with ones 
obtained via a different approach proposed by Narasimhan [24] and 
Ekelman [25]. The third example is the same example used by Williams 
[16]. Narasimhan modified the partial differential equation derived for 
an offset dual-reflector by Lee [19] to a ordinary differential equation 
which has to be solved numerically. 

6.2. Examples of Shaping of Three-Dimensional Dual-Reflector Antennas 
6.2.1. Example 3D1 

In this example, a shaped dual-reflector with a concave 

subreflector is to be generated so that a uniform aperture distribution 

is obtained at the aperture of the main reflector. The primary source 

has a cos^0 symmetric power pattern and is located 0.5' from the vertex 

of the shaped main reflector. The diameter of the main reflector is 

2.0'. The vertex of the shaped subreflector is located 0.2667' from the 

primary source. The maximum illumination angle of the primary source on 

the subreflector is 9 -31.42°. The initial points for the shaping are 

max 

thus given by 
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Figure 6.3. Shaped dual-reflector of example 3D1. 
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Figure 6.4. Shaped dual-reflector of example 3D1 obtained by the method 
of Narasimhan [24] and Ekelman [25]. 


and ( Z C0 »P C0 ) ^ determined from the ray condition given in Appendix C. 
The resulting shaped dual-reflector is shown in Figure 6.3. The 
surfaces obtained by using the method proposed by Narasimhan [24] and 
Ekelman [25] are shown in Figure 6.4. The results are in very good 
agreement as can be seen by comparing Figures 6.3 and 6.4. 

6.2.2. Example 3D2 

This example is similar to Example 301 except that the subreflector 

is a convex surface and the primary source is located 0.1' from the 

vertex of the main reflector. The maximum illumination angle 0 is 

° max 

18.26°, and the vertex of the subreflector is 0.4667' from the primary 
source. The initial points on the reflectors are 

(z^q, p^q) = (0 . 5667 ' , 0. 0) 

^ z 20’ p 20^ = 

and ( Z C0 >P C0 ) is thus determined. The resulting shaped dual-reflector 
is shown in Figure 6.5 and is again in very good agreement with the one 
obtained by using the method of Narasimhan [24] and Ekelman [25] which 
is shown in Figure 6.6. 

6.2.3. Example 3D3 

In this example, the conventional Cassegrain reflector used by 
Williams [16] is used as the starting point of shaping. The parameters 
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associated with this Cassegrain antenna are 


F =2.5' 
m 

F =1.5' 
c 

L =0.20943' 
v 

D =10.0' 
m 

84 5 

and a point source with cos * 0 power pattern is used as the primary 
feed. The secondary aperture is required to have uniform amplitude and 
phase distributions. The shaped surface obtained by Williams are shown 
in Figure 6.7. This is a enlarged version of Figure 3 of [16] . These 
surfaces are obtained by shaping the subreflector first and the main 
reflector is designed to create the desired uniform phase distribution. 
Figures 6.8 and 6.9 show the shaped reflectors calculated by the method 
of Narasimhan [24] and Ekelman [25]. 

In order to use the approach of shaping discussed in this chapter, 
the initial points for the shaping are chosen as 

(z 10 ,P 10 ) = (2. 18357', O') 


for the subreflector and main reflector, respectively. These points are 
actually obtained from Figure 6.7 which are the intersection points of 
the axis of symmetry and the shaped reflectors designed by Williams. 

The initial point of the caustic surface is thus determined. The 
resulting shaped surfaces are shown in Figures 6.10 and 6.11. These 
results show that very good agreement has been obtained among the three 
methods. 

6.3. Summary 

In this chapter, the shaping process for the two-dimensional dual- 
reflector was modified slightly for shaping of the three-dimensional 
dual-reflector with circular symmetry. Examples were used to verify the 
process. Comparisons were shown between different approaches. Although 
the method of synthesis discussed in this chapter is not much simpler 
than Narasimhan's approach, one can still examine the caustic surface of 
the subreflector reflected field as in the two-dimensional case. The 
behavior of the shaped reflectors can be explained physically. 
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Figure 6. 



■ 


. Shaped 3-dimensional circular symmetric dual-reflector from 
[16]. 
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Figure 6.8. Shaped dual-reflector of example 3D3 obtained by the method 
of Narasimhan [24] and Ekelman [25]. 
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Figure 6.9. Shaped subreflector of example 3D3 obtained by the method 
of Narasimhan [24] and Ekelman [25]. 
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Figure 6.11. Shaped subreflector of example 3D3 obtained by the new 
approach. 
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CHAPTER VII 


SUMMARY AND CONCLUSIONS 


A new approach for the shaping of dual-reflector antenna systems to 
generate a prescribed aperture distribution is proposed and studied in 
this research. The approach is based on the geometrical optics 
properties of conventional Cassegrain and Gregorian reflector antennas. 

A set of algebraic equations has been obtained and solved by a numerical 
method. This approach also explains how the shaped dual-reflector works 
for generating a prescribed aperture distribution by examining the 
behavior of the caustic curve of the shaped subreflector reflected 
fields. Several two-dimensional examples are presented to validate the 
formulation and solution of this new approach for shaping dual-reflector 
antennas. A conventional Cassegrain or Gregorian reflector antenna is 
provided as the starting point of shaping when a primary source power 
pattern is given. The method of physical optics is used to calculate 
the far-zone scattered fields from the subreflector and the near-zone 
scattered fields from the main reflector. The results are compared 
between the conventional and shaped reflectors. As expected , the shaped 
reflector system provides the aperture field prescribed. 

Rolled edge treatment to the shaped dual-reflector is also studied. 
The addition of the rolled edges to the shaped main reflector can reduce 
the edge diffracted fields and consequently, the ripple in the near-zone 
scattered fields. This improvement is especially important in the 
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application of compact range reflectors because creating a uniform plane 
wave illumination is a strict requirement. Corrections are made for the 
false end-point contributions of the rolled edges which are caused by 
the sharp termination of physical optics currents. By removing the 
false end-point contributions from the physical optics integration 
results, a more accurate representation of the scattered fields from 
reflector surfaces is obtained. Examples are also presented to show the 
improvement of adding blended rolled edges to the shaped main reflector. 

The shaping of three-dimensional circularly symmetric dual- 
reflectors is achieved by a simple modification of the two-dimensional 
case. The reflectors are Illuminated by primary sources with circularly 
symmetric power patterns and the secondary aperture distributions also 
have the property of circular symmetry. Examples are presented and 
compared with the results obtained by other approaches. 

Throughout this study, a numerical method is used to approximately 
solve for the surface equations. Different numerical methods could be 
used to solve these equations. 

Although some research has been done in the past for the shaping of 
dual-reflector antennas, only the shapes of the reflector surfaces were 
examined by most people. By using the method of shaping proposed in 
this study, one can clearly understand the detail of how the shaped 
reflector works to transform a given primary feed pattern into a 
prescribed aperture distribution. This is accomplished by examining the 
behavior of the caustic curve for the fields reflected from the shaped 
subreflector. Consequently, some potential problems associated with the 
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shaped reflectors can be disclosed through unusual properties of the 
caustic curve, such as a discontinuity in its slope. These possible 
problems affect the performance of the shaped reflectors especially in 
the compact range applications. 

Further study is recommended for applying the shaping approach 
proposed in this research to three-dimensional offset dual-reflectors. 
Unlike the circularly symmetric case for which the shaping can be done 
on a two-dimensional basis, the offset dual-reflector has to be shaped 
on a three-dimensional basis. It is pointed out in [14] and [26] that 
the exact solution to the dual-reflector synthesis does not exist in 
general. Consequently, an approximate solution must be sought when an 
exact solution cannot be found. Further investigation will have to be 
done on this topic. 
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APPENDIX A 


DIRECT APPLICATION OF PRINCIPLES OF GEOMETRICAL OPTICS 
TO DUAL-REFLECTOR ANTENNA SHAPING 


In this appendix, the principles of geometrical optics mentioned in 
Chapter III are used directly to generate a two-reflector antenna system 
from which a prescribed secondary aperture distribution is obtained. 

A.l. Formulation 

Consider the geometry shown in Figure A.l where a given primary 
source is located at the coordinate origin. Assume that the points 
<*10,yio> on the subreflector and on t ^ ie ma * n reflector are 

also known. For a given ray y=m^x which emanates from the primary 
source, the corresponding reflected ray from the main reflector, y=Y2> 
is to be generated and is parallel to the x-axis. The unknown point 
(X^,Y^) on the subreflector is assumed to be located on a line which 
passes through ( x jQ>y^o^ an< * ^ as t * ie s ^°P e Similarly, the unknown 

point ^2^2) on the main reflector is located on a line which passes 
through ( x 2 Q'y 2 (p an< * ^as s -t°P e m 2* These assumptions are valid when 
the surface points are very close to each other. From the application 
of the principles of geometrical optics and the above assumptions, a set 
of algebraic equations is generated and solved to obtain the surfaces. 
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Figure A.l. Geometry for direct application of principles of 
geometrical optics to dual-reflector shaping. 
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At the first reflector (i.e. the subreflector), which is given in 
Figure A. 2, the incident ray, reflected ray, and surface normal n^ 
satisfy Snell's law of reflection, in other words, 

e -6 = 0.-0 (A. 1) 


or 


20, =0+9 
1 o c 


(A. 2) 


By taking tangent on both sides of Equation (A. 2), one finds that 


tan(20j) = tan(0 Q +0 c ) 


(A. 3) 


From the identity 


o\ tana + tan|3 
tan< “ +8) * 1 tanatang 


Equation (A. 3) results in 


2tan0^ 

l-tan^0. 


tan0 +tan0 
o c_ 

l-tan0 tan0 
o c 


(A. 4) 


but , 
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PR IM ARY 
SOURCE 


Figure A. 2. Snell's lav on the subreflector. 


9 


tan© * m 
o o 


tan© = m 
c c 


tan©! 


1_ 

m l 


Consequently, Equation (A. 4) becomes 


-2m^ 

m ^-1 


m +m 
jo 

1-m : 


c 

m 


o c 


(A. 5) 


For the main reflector, as shown in Figure A. 3, Snell's law states that 


r 


e 

c 


20„ 


in which 


e 


c 



(A. 6) 


(A. 7) 


Substituting (A. 7) into (A. 6) and taking the tangent on both sides of 
the equation, the following equation is obtained: 
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Figure A. 



to' 


3. Snell's law on the main reflector. 
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Also notice that X^, have to satisfy 


Yj - m Q X 1 (A. 9) 

Y l _y 10 “ m l( X r x 10> (A. 10) 


and X2 has to satisfy 

Y 2 _y 20 = m 2^ X 2 _X 20^ (A. 11) 

In the above formulations, m Q and Y2 are assumed known so there are 
six unknowns Xj,Y^,X2,mj,m2,m c , with six equations which are summarized 
as follows: 


- 2 m^ 

TT 

m^-l 


m +m 
o c 

1-m m 
o c 


m 


c 


- 2 m 2 

TT 

m 2 1 


m 


c 



(A. 12) 


(A. 13) 


(A. 14) 


Y- = m X- 
1 o 1 


(A. 15) 


Y i y io _ m i< x r x io> 


(A. 16) 
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It is apparent that the above equations have to be solved numerically in 
order to obtain the surface of the reflectors. In the next section, a 
perturbation method is presented to solve the above equations. 


A. 2. Method of Solution 

From Equations (A. 15) and (A. 16), it is found that 


X 


1 


y io Vio 

m -m- 
o 1 


(A. 18) 


Y 


1 


m 

<y io-"i x io> 

O 1 


(A. 19) 


From Equation (A. 17), 

X 2 " jq < V Jr 20 + V20 ) < A ' 20 > 

Substituting (A. 18), (A. 19), and (A. 20) into (A. 14), the following 
equation is obtained, 



Equations (A. 12), (A. 13), and (A. 21) involve the unknowns m.,m 0 , and m 

1 l c 

and are very difficult to solve exactly. Consequently, a perturbation 
method which assumes that 


m 


c 


m 


co 


+ Am 

c 


(A. 22) 


m l = m 10 + ^"l (A. 23) 

m 2 = m 20 + ^ m 2 (A. 24) 

is used to solve m^, m 2 , and m £ numerically. In Equations (A. 22) 
through (A. 24), m co , m^, and ii^q are the corresponding slopes m c> m^ 
and m 2 for the previous ray and surfaces; Am c , Am^, and Am 2 are small 
perturbation of m £ , n^, and n^. By substituting these three equations 
into Equations (A. 12), (A. 13), and (A. 21) and neglecting the higher 
order terms such as An^Am^, An^An^, Am^Anu, etc., the following three 
linear equations in Am c , Am^, and An^ are obtained. 


A ll^ m c + A 12 Am l + A 13^ m 2 " B 1 


A 21*"c + ^2^1 + A 23 Am 2 = B 2 
A 31 to c + A 32 Am l + A 33 Am 2 = B 3 


(A. 25) 
(A. 26) 
(A. 27) 


where 
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A 11 “ 1 + 2m o m 10 " m 10 
A 12 = 2 ( m o m co -m o m 10 -m 10 m co -1 ) 

A 13 = ° 

B 1 * 2m io( 1 - m 0 m co) + (V m co) ( m 10- 1 ) 

a 2 , 

A 21 * m 20 _1 

A 22 = 0 

A 23 " 2 ( ra co m 20 +1 ) 

B 2 “ m co( 1_m 20) ~ 2m 20 

A 31 = ( m o -m lo) ( Y 2 _y 20 +m 20 x 2o) “ m 20 ( y 10 -ra 10 X lo) 

A 32 " m co[ ra 20 x 10 - (V y 20 +m 20 x 2o)] + m 20 (Wio) 

A 33 " m co[ X 20( m o- m 10)-( y 10- m 10 X 10)]-[( m o- ra 10) Y 2 " m o ( y l(T m 10 x lo)] 
B 3 = [ m 20(V m 10) Y 2 ' m o m 20 ( y 10 -m 10 X lo) ] 

" “co [ ( m o -m lo) ( Y 2 -y 20 +m 20 x 2o) ~ m 20 ( y 10~ m 10 x lo)] 


185 



The above equations are solved for Am c> , and Am 2 as 


Am 

c 


A 12 A 23 B 3 A 23 A 32 B 1~ A 12 A 33 B 2 
A 31 A 12 A 23 _A 11 A 23 A 32 _A 12 A 21 A 33 


(A. 28) 


- I^ 2 ( B l- A ll to c) (A- 29 ) 

*"2 ■ iq; ( B 2- A 21 to c) < A - 30 > 


and the unknowns x i> Y i> X 2 » m i» m 2 » m c are obtained. 

It is noted that the above method is valid provided that m^, n^Q, 
n» co are known and Am^, Anig, Am^ are very small. There is a difficulty 
associated with this approach. When the initial slopes are very large, 
the changes of the slopes Am^, An^, or Am c are not small and the 
solutions are not accurate. Another method for solving the surface 
equations might have to be used. 


APPENDIX B 


NUMERICAL METHOD FOR SOLVING THE SURFACE EQUATIONS 

As discussed in Chapter III the equations related to the shaped 
reflector, i.e., Equations (3.1) through (3.6) in that chapter, can be 
reduced to 

Vc- y c[(v s, c) 2t4p v 4p2 ]; 

- IWcKMf' *cV y c[(V y c) 2 * 4F x c- 4r2 ] ; • 

( 4F (VWcK [( Y 2- y c) 2 * 4F V 4f2 ] ’ 

- - 5 K- 4 ,! ) ! { w (vwJ- , .[W 2,4p v 4p2 ] • 

(B.l) 


where all the variables have been defined in Chapter III. In Equation 

(B.l), m and Y 0 are known and, F and F are functions of x and y 
o z c c c 

only. In addition, x and y are related to x and y by 

c ■'c co ' co ' 


» 


x = x + cosa Ap 
c co o 


= y + sina Ap 
J co o 


(B.2) 
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where Ap is the unknown distance between (x ,y ) and (x ,y ) and a is 

CO CO c c o 

the angle of the old ray with respect to the positive x-axis which is 
defined in Chapter III. 

By substituting Equation (B.2) into (B.l), an 8 th order polynomial 
equation in Ap has to be solved numerically. Although one can try to 
solve Ap numerically directly from the polynomial equation, it is still 
very tedious. A method which assumes Ap is small is used in this study 
to solve Ap. This assumption is an appropriate one since that the 
incremental ray tube formed by the old and new rays is very small so 

that one will expect that variation of Ap is very small. Consequently, 

2 3 

the higher order terms such as (Ap) ,(Ap) , ...,etc., in the polynomial 
equation are neglected so that a linear equation in Ap is obtained and 
solved. In the following, different terms in Equation (B.l) are 
discussed separately by substituting Equation (B.2) for x £ and y c into 
Equation (B.l). Also notice that 

P . \ (l-sec« 0 ) (x c -x 2Q ) 

or 

F - - 1 tan r (v y 2o) 


The first expression of F is used in the following analyses. 
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(1) ( Y 2- y c) 2 + 4F X C - 4p2 

- ( y co + sinot o Ap - y 2 ) 2 + 2 ( 1 - seca o)( x co + cosa o Ap - x 2o) 

(x co + cosc^Ap) - (l-seca o ) 2 (x co + cosc^Ap - x 20 ) 2 

= ( y co" Y 2) 2 + 2 sina o( y co" Y 2) Ap + sin 2 a Q (Ap) 2 + 2 (l-secaj • 

[ X co ( X co -X 2o) + ( X co- X 2o) cosa o Ap + X co cosa o Ap + cos2 « 0 < A P> 2 ] 

- [ 1-seca 0 ] 2 [( x C o _X 2o) 2+2coSa o ( x co -X 2o) Ap + ^os 2 a o (Ap) 2 ] 

" {( y co- Y 2) 2 + 2 ( 1_seca o) x C o ( x co _x 2o) - {'-^“oYKo-'lo) 2 } 

+ 2 { sina o( y co- Y 2) + ( 1-seca 0 ) ( 2x co -x 2o) COSO o - cos * o [ 1 - seca o ) 2 * 
( x co -X 2o)} Ap + { sin2 V 2 cos2 «o( 1 - sec %) - cos 2 a o (l-seca o ) 2 }(Ap)‘ 


It can be shown that the coefficient of (Ap) is zero. Thus, 


(V y c) 2 + 4F V 4p2 " Z 1 + Z 2 Ap 


(B.3) 


in which 


h ’ (*co- Y 2) 2 * 2 (l-“'«o)*co ( x co- x 2o) - ( 1 - sec “o) 2 ( x co- x 2o)' 
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Z 2 . 2{sin« o (y co -Y 2 ) * (cos« 0 -l) ((l^ecajx^-x^sec.j 


(2) 4a 


- (&♦ y 10 ) 1/2 ± [( x 10- x c) 2 * (^lO-yc) 2 


2 +: ellipse 

’ hyperbola 


ft 


x 2 + y 2 ] 1/2 + 

X 10 + y 10j - 


x ln -x 
10 c 


1 + 


frl0- y c^ 


l x 10 -x cj 


2 ' 1 / 2|2 


but 


y 10‘ y c , 

X -x = tan °o 

x 10 c 0 


Thus, 


4a 


fx 2 + y 2 ) 1/2 + 

l x 10 + y ioJ * 


x,„-x 
10 c 


seca 


when seca Q <0, 


x^q-x c <0 for hyperbola 


x^q-x >0 for ellipse 


when seca Q >0, 
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for hyperbola 


x^q-x c <0 for ellipse 
Consequently, 


4a2 ■ { ( x 10 + y lo) 1/2+ ( x c _x lo) seca o} 2 


The geometry for the above analysis is shown in Figure B.l. 


2 2 
(3) F* - 4a Z 


(2 21 f 

( 2 2 U/2 | 

\ 

f ) 

l x c * y cJ - [ 

L x io + y ioJ * 

x -x. n seca > 
l c 10; o 

/ 

[( x co-' coso o“ > ) 

I 2 * ( y co* si "V p ) 

I 2 ] 


[(*10 * 4) 1/2 + ( x co~ x lo) seca o] * cos V ec “o 4p 


2 2 2 2 
x + 2 x cosa Ap + cos a (Ap) + y + 2 y sina Ap 
co co o o co co o 


sin 2 « o (4p) 2 - [(* 2 0 * y 2 0 ) 1/2 + (x co -x 10 ) S .c« o ] 


- 2 [( x 10 * y lo) 1/2 * ( x co- X lo) Sec vb - 
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Z 3 + Z^Ap 


in which 


CM 

CM 

7 2 2 U/2 

( ) 12 

= x +y - 
co •'co 

O 

+ 

<< 

o 

+ 

l x co“ x loj seco oj 


• 2 { x co cos V y co sin v[( x 10 * y lo) 1/2t ( x co- x lo) sec “o]} 


(4) x +m y = [x +cosa Apl + m [y +sina Apl 
c o' c k co o ) o v co o ) 


— + Z 6 Ap 


where Z* = x + m y 
5 co 


co 


Z. = cosa + m sina 
6 o o o 


(5) Prom (2), 4a 2 . {(*i 0 +yi 0 ) 1/2 * (v x lo) secot o} 2 
‘ {[( x 10* y lo) 1/2 + ( x co -x lo) secot o] * ‘“’F 
■ [( x 10 +y lo) 1/2 ♦ ( x co -x lo) secot o] 2 
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+ 2 [( x 10* y lo) 1/2 + ( x co -x lo) seca o] + < 4 ’> 
Z 7 + Zg Ap 

where 

Z 7 - [( x 10* y lo) 1/2 + ( x co~ x lo) secot o] 2 
Z B * 2 [( x 10 +y lo) 1/2 * ( x co -x lo) secot o] 

(6) 4Px c Y 2 . 2 (l-secaj (v*2o) x c Y 2 

‘ 2 ( 2 ~ seca o) y 2 ( x co + C0s “o 4 »- X 2o) ( x co +cosa o‘ 
' 2(1 - seca o )Y 2{ x co( x co- x 2o) * ( 2 x c<T x 2o) co! 

+ cos^a^Ap)^ 

~ z 9 + Z 1Q Ap 

where 


Z Q = 2 fl-seca |Y- x [x -x„ n l 
9 V oj 2 col co 20) 



Z l0 - 2 (l-sec« 0 )Y 2 (2x co -x 20 ) 


cosa 


(7) 4F(Y 2 -y c * m o xJ 


■ 2 W ( x c -x 2o) [V y c * Vt] 

■ 2 (‘-"“o) [( x co~ x 2o) + cos “o ap ] ( Y 2*Vco- y co) 
+ |m o cosa o -sina^Apj 


= 2 


+ 


(i-sec« 0 ) { ( x co - x 2o) (V Voodoo) 

[( x co” x 2o) (v os “o “ s in “o) » cos “o ( Y 2 *"o x co- y co) ] 4p 


+ cosa 


m cosa -sina 
4 o o o 



_ Z ll + Z l2 Ap 


where 


Z 11 - 2 ( l -“ c «o) ( x co- x 2o) (Wc<>- y co) 


Z 12 * 2 ( 1 - seca o)[( x co- x 2o)(v os “o- sln “o) tcosa o( Y 2*Vco- y co 
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^x c+ m o y c j -^a 

■ (W ’) 2 - (W>) ( U "o) 

' [ Z 5 - Z 7 ( 1+B o)] + [ 2Z 5 Z 6 - KK] 
= W 1 + W 2 Ap 

where 

"l - Z 5 - Z 7 ( U "o) 

h , 2Z 5 Z 6 - (l« 2 )z e 

(9) 4FX C Y 2 - y c [( Y 2 -y c ) 2 ♦ 4FX C -4P 2 ] 

- V Z 10 6e - ( y co tsin “o 4p ) [ Z l* Z 2 4p ] 

- Z 9 +Z 10^ p - [ y coV ( y coV p *"«o Z i) 4 f 
' [V y co Z l] + { Z 10 - ( y co Z 2 + sln “o Z l 


+ sina Z 
o 

)h 


= V. + W.Ap 



in which 


W 3 - V F c 0 2 l 

“4 ■ Z 10 - y co Z 2 - sin “o Z l 

do) 4F(v 2 -y c *» 0 x c ) - " 0 [(vy c ) 2 * 4F,I c - 4p2 ] 
* Z 11 * Z 12 4p - "o( Z l + Z 2 4 ') 

" ( Z U ' Vl) * ( Z 12'"o Z 2) ip 
■ *5 * V 

in which 

W 5 " Z 11 - m o Z l 
W 6 = Z 12 " m o Z 2 

By substituting (1) - (10) into Equation (B.l): 
1 st terra of lefthand side of (B.l): 
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" l V l + V P J ( W 3 + V P J 


= (w i + W 2 Ap] [W3 + 2W 3 W h L 

* V 1 + V p 

e 

V 1 - B 1 U 3 

V 2 . 2W 1 »3» 4 * v 2 «2 

term of lefthand side of (: 

- (W") (VV*) (V B , 

' VV P 

ft 

V 3 = Z 3 Z 5 W 3 W 5 

= z,z« (v,V A+ V y V J + W 0 W 


+ V*(Ap) 


'• 1) : 

Ap ) ( W 5 + V 


c fz_z,+z,z. 


Righthand side of (B.l): 


- - i (W’WVV^ 


W’ 


where 


V c = - 


1 2 2 
4 Z ?5 


V 6 - - I ( Z 3 Z 4 W 5 + Z 3 U 5 U 6) 


Finally, from Equation (B.l) 


(W») - (W<“) * VV’ 


thus 


Ap = 


■Wsi 

VV V 6 


From Ap, the surfaces of the reflector and the caustic curve are solved. 
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APPENDIX C 


GEOMETRICAL OPTICS REFLECTED FIELDS 
FOR CASSEGRAIN AND GREGORIAN REFLECTOR ANTENNAS 


In this appendix, the geometrical optics reflected fields for the 
Cassegrain and Gregorian reflector antennas are calculated. In 
addition, the conditions which the initial points on the subreflector, 
main reflector, and caustics curve for the shaped dual-reflector must 
satisfy are derived in this appendix. 


C.l. Reflected Fields for the Two-Dimensional Cassegrain Reflector 
Antenna 

The geometry of a Cassegrain reflector antenna is given in Figure 
2.6 of Chapter II. As discussed in Chapter II, for a source located at 
the real focus of the subreflector, the reflected field from the 
subreflector appears to emanate from the virtual focus. In other words, 
the virtual focus is the caustic of the subreflector reflected fields. 
Assuming that the primary source is a magnetic line source, the 
reflected magnetic field is calculated as 


H Z ( Sl ) = H^) 





1 


-jkSj 

e 


(C.l) 
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in which H^(Q^) is the incident magnetic field at the point of 

reflection on the subreflector, p£ is the reflected field caustic 

distance and s^ is the distance from the point of reflection to the 

£ 

field point. In this case, is the distance between the virtual focus 
and the point of reflection, i.e., Pj=P c as shown in Figure C.l. The 
far-zone reflected field is calculated by letting s^-»“ and is given as 


Hg(s-*») = HjCQj) 



(C.2) 


The reflected field from the parabolic main reflector is calculated by 


H m< s 2> 


Hj(Q 2 ) 


r 

p 2 


-jits. 


(C.3) 


"2 * s : 


in which H*(Q 2 ) is the incident magnetic field at the reflection point 
Q 2 , p 2 is the reflected field caustic distance, and s 2 is the distance 
from the point of reflection to the point of observation. For a 
focussed parabola, the reflected field caustic distance p 2 is infinite 
and the reflected ray is parallel to the axis of symmetry. 
Consequently, 


H >2> " "iCV^ 2 


(C.4) 
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Figure C. 


OBSERVATION 



. Geometrical optics reflected field from a hyperboli 


subreflector. 


For the Cassegrain reflector, the incident field H?(Q 2 ) is the reflected 
field of the subreflector and is given by 


Hj(Q 2 ) = H^Qj) 


. P +s i 

sc 1 


-jks 1 

e 


(C.5) 


where s^ is the distance between and Q 2 as shown in Figure C.2. 
Thus, the reflected field of the parabolic reflector is given by 


H m< s 2> « H i<°l> 


P c -jk(s 1+ s 2 ) 


p c +s l 


(C.6) 


If the magnetic line source has a far-zone pattern of F($), 


H^Qj) = |F(4^) 


i e 


-jkp. 




(C.7) 


Thus, 


“> 2 > ’ F*? 


-jkp. 


FT 


P c -jk(s 1 +s 2 ) 


P c +S l 


(C.8) 
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y 



Figure C.2. Ray geometry of a Cassegrain reflector for calculation of 
reflected fields. 


204 


C.2. Reflected Fields for the Two-Dimensional Gregorian Reflector 
Antenna 

The geometry of a Gregorian reflector antenna is given in Figure 
2.8 of Chapter II. For a line source located at the first focus of the 
elliptic subreflector, the reflected field from the subreflector appears 
to emanate from the second focus. Thus, the second focus is the caustic 
of the subreflector reflected field. The reflected field from the 
subreflector is calculated by Equation (C.l), except that the reflected 
field caustic distance is negative since the caustic is located in 
between the points of reflection and observation as shown in Figure C.3. 
Thus, 


H s< s l> 


H i«V 


-P„ 


V s i 


-jks x 

e 


(C.9) 


where p c is the distance between and the caustic and is positive. 
The far-zone reflected field is given as 



(C.10) 


Similar to the Cassegrain case, the reflected field from the parabolic 
reflector is given by 
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Figure C.3. Geometrical optics reflected field from an elliptic 
subreflector. 
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(C.ll) 


-jkp. 


H m< S 2> - JF(V 


n e 




-p„ 


- p c* S 1 


-jk(s 1 +s 2 ) 

e 


The ray geometry of this Gregorian antenna is shown in Figure C.4 where 

s, is the distance between Q. and Q„ and p is the distance between Q 1 
1 1 l c 1 

and the second focus. 


C.3. Relationship of Initial Points For Shaping of Two-Dimensional 
Dual-Reflector Antenna 

It is mentioned in Chapters III and IV that in order to solve the 
surface equations, a set of initial points on the subreflector, main 
reflector, and caustic curve must be provided a priori. For given 
primary source and secondary aperture distributions, these initial 
points have to satisfy a ray condition. This condition is derived in 
this section for a dual-reflector with either a convex or concave 
subreflector. 

C.3.1. Dual-Reflector with Convex Subreflector 

The basic assumption in this research for dual-reflector shaping is 
that the corresponding points on the subreflector, main reflector and 
caustic curve are points on a Cassegrain or Gregorian antenna. 

Connecting these points forms a ray such that the reflected field at the 
corresponding aperture point is given by 
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Figure C.4. Ray geometry of Gregorian reflector for calculation of 
reflected fields. 
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(C. 12) 


h z = mr 

m v ' 


-jkp. 


J'T 


V p c 


-jk(s 1 +s 2 ) 

e 


when the subreflector is a convex surface. Assuming that the aperture 
distribution corresponding to this ray is known as U q , then 


-jkp. 


UFOfr) 


i e 


-jk(s 1 +s 2 ) 


S l +P c 


= U 


(C. 13) 


which simplifies to 


F(<fr) 


1_ 

Pa 


= U 

P c + Sj O 


(C. 14) 


This is the condition that each ray has to follow. In Equation (C.14), 
pj, p c and Sj^ are functions of (Xj.y^, (x 2 ,y 2 ) and (* c> y c ) and U Q is a 
function of y 2 - The initial points ( x 10 ' y 10^’ * X 20’ y 20^ and ^ X co ,y co^ 
have to follow Equation (C.14) also. In this study, it is assumed that 
the shaping process starts along the primary feed axis which is tilted 
by Consequently, one can assume that | F( *o ) I = 1 and we obtain 


1_ _JiL_ 
p i p c +s l 


U o (y 20> 


(C.15) 
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Although there are many ways to determine the initial points from 
(C. 15) , it is noted that y 2Q has to be known in order to determine U q 
and also ( x io ,y 10^ has t0 ^ e on the P r i mar y source axis. Thus, 

< x 10 ,y 10> and (x 20’ y 20^ are assumed and ^ x co ,y co^ * S s0 ^ ved f rom (C.15), 
as shown below. 

Since ( x jo ,y 10^ and ^ x 20’ y 20^ are known ' p i and s i are known for 
the initial ray. From Equation (C.15) one finds that 



(C. 16) 


From the geometry shown in Figure C.5, it is found that 


X CQ- X 10 y CQ- y iO ^C 
X 10 -X 20 ' y 10- y 20 ' S 1 


(C. 17) 



From these initial points, one can start shaping the reflector above and 
below the ray which is formed by connecting these initial points. 
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Figure C.5. Relation of initial points for a shaped dual reflector with 
a convex subreflector. 

y 



Figure C.6. Relation of initial points for a shaped dual reflector with 
a concave subreflector. 
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C.3.2. Dual-Reflector with Concave Subreflector 
Similar to the convex subreflector case, the corresponding points 
on the subreflector, main reflector and caustic curve for this case are 
on a Gregorian reflector. The reflected field at the corresponding 
aperture point is given by 


H z = -J F( «|>) ' 
m 




s r p c 


-jk(s 1+ s 2 ) 

e 


where p >0. Thus, one finds that 
c 



-P -ikp. 

K c J i 
e 

Vc 


-jk(s 1 +s 2 ) 

e 


2 



(C. 20) 


(C.21) 


or 




(C. 22) 


Again, by assuming that the shaping starts along the feed axis with a 
tilt angle 4> 0 and |F(+ o )|=l, the following relationship is obtained: 


1_ 

Pi 


s l _p c 



(C. 23) 
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From the same assumption that ( x ^o ,y 10^ an< * ^ x 20 ,y 2(P are so t * iat 


and Sj are known, then 


P 

s 


c 

1 


p i U c 


1+ Pl U o 


(C. 24) 


and 


X CQ- X 10 y CQ- y iO 
x 20~ x 10 " y 20 -y 10 ' S 1 


(C.25) 


Thus, 


co 


p i U o 

1+p i U o 


< x 10' x 20 ) + x 10 


(C.26) 


p i U o 

y co = - UTIT (y 10- y 20> + y 10 
i o 


(C.27) 


C.4. Relationship of Initial Points for the Shaping of A Three- 
Dimensional Circularly Symmetric Dual-Reflector 
The major difference between the two-dimensional and three- 
dimensional case is that the spread factor of the reflected field in the 
three-dimensional case is dependent on the two radii of curvature since 
the ray tube is three-dimensional. By simple modification of the two- 
dimensional case, the condition which the initial points have to satisfy 
for the three-dimensional case can be obtained. 
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C.4.1. Dual-Reflector with Convex Surface 

For a point source with radiation power density F(6,<f>) illuminating 
the subreflector of the dual-reflector, the reflected field at the 
aperture of the main reflector is calculated by 


u r = >lF(e,<fr) 


1 e 


-jkp. 


P, P 


12 


-jk(s 1 +s 2 ) 


(p 1 +s 1 )(p 2 +s 1 ) 


(C.28) 


in which p. is the distance between the primary source and the point of 
reflection on the subreflector, pj and p 2 are the subreflector reflected 
field caustic distances, s^ is the distance between the point of 
reflection of the subreflector and main reflector, and s 2 is the 
distance from the point of reflection on the main reflector to the 
aperture. 

For the shaped dual-reflector with a convex subreflector, the 
corresponding points on the subreflector and main reflector are assumed 
to be located on a particular Cassegrain reflector; consequently, the 
two radii of curvature p^ and p 2 are the same and are equal to the 
distance p £ between the caustic and point of reflection on the 
subreflector, thus, Equation (C.28) becomes: 
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r 


-jkPj 

4f(m5T !L 7 — 


l p c +s U 


-jk(s 1 +s 2 ) 
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(C.29) 
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It is assumed that the power density at the aperture point is given by 


U , and is related to u as 
o r 


-jkpj 

^F<e,$) ‘ 

p cl 

-jk(s 1 +s 2 ) 

e 

p i 

[ p c +s lj 



(C.30) 


Assume that the shaping begins from the axis of symmetry, thus, let 
F(0,$) = F(0,0)=1, one finds that 


1 _ 

Pa 


V S 1 


■K 


(C.31) 


Comparing Equation (C.31) with (C.15), the condition that the initial 
caustic point has to satisfy is given by 




(p 10 p 20 ) + p 10 


(C. 32) 
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CO 



(z 10 -z 20 ) + z 10 


(C.33) 


provided that ( z jo ,p 10^ anc * ^ z 20 ,p 20^’ initial points on the 
subreflector and main reflector are known. Note that the coordinates in 
the three-dimensional case are (z,p). 
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C.4.2. Dual-Reflector with Concave Subreflector 
The difference between this case and the convex subreflector case 
is that the subreflector reflected field caustic distance is negative, 
i.e. , 


P 2 - -P c 


(C.34) 


in which p £ is the distance between the caustic and the point of 
reflection on the concave subreflector. Thus, from Equation (C.28), it 
is found that 
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-jkpj 
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(-Pc+S!) 


-jk(s 1 +s 2 ) 
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(C.35) 


Consequently, 


-jkpj 

mrw S -r— 
p i 


-P, 


" p c +s l 
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or, by letting F(0,<|>) = F(0,0)=1, 


(C.36) 


1 



(C. 37) 


Thus, 
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